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We derive the spin-atomic vibration interaction Vsa and the spin-flip Hamiltonian Vsf of 
a single atomic spin in a crystal field. We here apply the perturbation theory to a model with 
the spin-orbit interaction and the kinetic and potential energies of electrons. The model also 
takes into account the difference in vibration displacement between an effective nucleus and 
electrons, Ar. Examining the coefficients of Vsa and Vsf, we first show that Vsa appears for 
Ar^O, while Vsf is present independently of Ar. As an application, we next obtain Vsa and 
Vsf of an Fe ion in a crystal held of tetragonal symmetry. It is found that the magnitudes of 
the coefficients of Vsa can be larger than those of the conventional spin-phonon interaction 
depending on vibration frequency. In addition, transition probabilities per unit time due to 
Vsa and V"sf are investigated for the Fe ion with an anisotropy energy of —\D\ S\ , where D 
is an anisotropy constant and Sz is the Z component of a spin operator. 

KEYWORDS: spin-atomic vibration interaction, spin-flip Hamiltonian, spin-phonon interaction, 
single atomic spin, spin-orbit interaction, crystal field, surface adsorbate, per- 
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1. Introduction 

Recently, the magnetic properties of molecular spin 1-9 and atomic spin 10-19 systems have 
been extensively studied to develop ultimate microscopic elements for mass-storage devices 
and quantum information devices. In particular, the spin relaxation 1-3 ' 19 has attracted much 
attention from the viewpoints of data writing and storage. This spin relaxation 20 means that 
the spin relaxes from the excited state to the ground spin state via interactions between 
the spin and other degrees of freedom. The degrees of freedom are considered to be, for 
example, the atomic and lattice vibrations. Here, the vibrational states reflect structures of 
the systems. Typical structures of the above atomic spin systems are a magnetic ion adsorbed 
on a substrate 10 or the magnetic ion trapped between a tip and a substrate. 21 ' 22 For such 
atomic spin systems, however, theoretical studies of the interactions between the spin and 
vibrations have scarcely been performed so far. 

A relevant interaction is considered to be spin-phonon interactions, 23-30 in which a local- 
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ized spin interacts with the phonon of a lattice system. The expressions of the spin-phonon 
interaction 23-29 have been originally derived using the spin-orbit interaction and the modula- 
tion of crystal field potential energy due to the lattice vibration. In particular, Van Vleck 23 ' 24 
obtained an expression applying a perturbation theory to the following model: the crystal field 
potential energy in an equilibrium state was included in the unperturbed Hamiltonian, while 
the spin-orbit interaction and the modulation of crystal field potential energy were the per- 
turbed one. 29 Mattuck and Strandberg 25 also derived essentially the same expression as Van 
Vleck's one 23 using a slightly different method. On the other hand, Hartmann-Boutron et a/. 30 
proposed a spin-phonon interaction applying the second-order perturbation theory to a model 
with the spin-orbit interaction. Here, the coefficient of the spin-phonon interaction contained 
the local strain on the spin site due to the phonon. 31 In these spin-phonon interactions, the 
structure of the systems was assumed to be a periodic structure (e.g., a bulk crystal), where 
a unit cell consisted of the magnetic ion and the surrounding ions. 32 The phonon then had a 
continuous energy spectrum on the assumption that the phonon frequency was proportional 
to the magnitude of the wave vector. 

Such spin-phonon interactions, however, may not be suitable for the atomic spin on sub- 
strate (or surface), because the vibration of the adsorbate is not described in their interactions. 
The vibration of the adsorbate on surface is obviously different from the above-mentioned 
phonon. The local vibrational density of states at the adsorbate shows sharp peaks at certain 
energies; 33 that is, this adsorbate has discrete vibrational energy levels. In contrast, the above 
phonon exhibits the broad vibrational density of states, 33 i.e., continuous vibrational energy 
levels. 

Furthermore, in the spin-phonon interactions, the difference in vibration displacement 
between an effective nucleus and electrons has not been taken into account. Here, this nucleus 
consists of a nucleus and core electrons. As shown in Appendix , the difference in displacement 
can be actually obtained, when the nucleus (plus charge) and electrons (minus charge) have the 
displacement in the presence of a crystalline electric field due to the surrounding ions (minus 
charge) (see Fig. 1). The displacements of the nucleus and the electrons are written by Ar n 
and (1 — rj)Ar n , respectively, where rj (0 < rj < 1) is a dimensionless quantity characterizing 
the difference. For instance, this r\ is evaluated to be about 0.05 for a single Fe ion on the 
CuN surface (see Appendix ). 

On the basis of the above-mentioned facts, we aim to derive a spin-atomic vibration 
interaction of a single atomic spin surrounded by several ions, where the ions are present in 
the substrate. Here, the spin-atomic vibration interaction represents an interaction between 
the spin of the single magnetic ion and the vibration of this ion. This interaction also takes 
into account the above rj. 

In this paper, we derived the spin-atomic vibration interaction and a spin-flip Hamiltonian 
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of a single atomic spin in a crystal field, applying a perturbation theory to a model with 
the spin-orbit interaction and the kinetic and potential energies of electrons. The spin-flip 
Hamiltonian contains spin- flip operators such as SxSz and SySz, but no atomic vibration 
operators, where Sj is the I component of the spin S=(Sx, Sy, Sz)- From the coefficients 
of the spin-atomic vibration interaction and the spin-flip Hamiltonian, we found that the 
former appears for ?77^0, while the latter is present independently of n. We also obtained their 
coefficients in the case of the Fe ion in a crystal field of tetragonal symmetry. The magnitude 
of the coefficient of the spin-atomic vibration interaction was estimated to be larger than that 
of the conventional spin-phonon interaction in a specific region of the vibration frequency. In 
addition, we discussed transition probabilities per unit time due to the spin-atomic vibration 
interaction and spin-flip Hamiltonian of the Fe ion. 

The present paper is organized as follows: In §2, we propose a model with the spin-orbit 
interaction and the kinetic and potential energies of electrons. This model also takes into ac- 
count the vibration displacement. In §3, we derive the spin-atomic vibration interaction and 
spin-flip Hamiltonian. In §4, we investigate the vibration frequency dependences of their coef- 
ficients of the Fe ion in the crystal field of tetragonal symmetry. In addition, their coefficients 
are evaluated by assigning appropriate values to the spin-orbit coupling constant. In §5, we 
compare the spin-atomic vibration interaction with the conventional spin-phonon interactions. 
In §6, we discuss the transition probabilities per unit time due to the spin-atomic vibration 
interaction and the spin-flip Hamiltonian. The conclusion is presented in §7. In Appendix , 
we evaluate the difference in vibration displacement between the effective nucleus and the 
electron. The vibration displacement of the nucleus is evaluated in Appendix . In Appendix 
, we describe the expressions of the coefficients of the spin-atomic vibration interaction and 
spin-flip Hamiltonian of the Fe ion. 

2. Model 

Toward the calculation of perturbation energies in §3 (i.e., the spin-atomic vibration inter- 
action and the spin-flip Hamiltonian), we here derive the unperturbed and perturbed Hamil- 
tonians for a model with the spin-orbit interaction and kinetic and potential energies of 
electrons. This model also takes into account the difference in vibration displacement between 
the effective nucleus and electrons. In addition, we introduce orbital and spin states for the 
unperturbed Hamiltonian. 

2.1 Hamiltonian in an equilibrium state 

We first propose a Hamiltonian in an equilibrium state (with no vibration displacement) for 
3d electrons of an iron-group ion (i.e., a magnetic ion) in a crystal field due to the surrounding 
ions. Here, the 3d shell of the magnetic ion is incompletely filled, and the surrounding ions have 
a minus charge. In particular, we treat the case of a weak crystal field. Namely, the electrons 
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Fig. 1. Positions and displacements of the effective nucleus (the large circle) and the £th electron 
(the small circle) in the single magnetic ion. The magnetic ion is surrounded by several ions. Here, 
one of the surrounding ions is shown by the medium circle, rf. position vector of the £th electron 
measured from the nucleus at an equilibrium position. Ar n : displacement of the nucleus. Ar^: 
displacement of the £th electron, r'f position vector of the displaced £th electron measured from 
the displaced nucleus. 

occupy energy levels from the bottom keeping their spins parallel, according to Hund's first 
rule. 34 The system is here assumed to have the total spin S. 

We now focus on electrons existing in the partially filled cr-spin shell of the magnetic ion, 
with a being up or down. When the number of electrons of the ground state is less than half 
(more than half), the cr-spin shell corresponds to the up-spin shell (down-spin shell). The 
excited states are then assumed to consist of orbitals in the cr-spin shell. In addition, the 
electrons in the cr-spin shell have the spin-orbit interaction V so . S5 Here, |V^ | is considered to 
be smaller than the magnitude of the crystal field potential energy. 25 When the number of 
electrons in the cr-spin shell is N, we have the following Hamiltonian for the electrons: 



n = n e + v s 



SO) 



(1) 





V so = XL-S, 



(3) 



with 



Vt = {pe,x, Pe,y,Pe,z) = mr e , 



(4) 




(5) 




(6) 




(7) 
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1 



N 



L — (L x , L y , L z ) 



h 



(9) 



withr £ =(a^, y e , z e ), r e =\r e \ =\jx\ + yj + zj, S=(S X , S y , S z ), C±= ±{^/^)g^ i Z cS /2S, and 
fi B =eh/2m. Here, % e denotes the Hamiltonian consisting of the kinetic energy J2eLiPj/(^ m ) 
and the potential energy ^2^ = i[V en (r£) + V c (r£)]. The quantity m is the electron mass, T(_ is the 
position of the ith. electron measured from the nucleus (see Fig. 1), and p e is the momentum 
of the ith. electron. The operator V en (ri) is a spherically symmetric potential energy of the 
ith electron created by a nucleus and the other electrons in the magnetic ion. The quantity 
— e (< 0) is the electronic charge, and eo is the electric constant. In addition, Z c g represents 
the effective nuclear charge for the 3d electron, 36 which can be evaluated from Slater's rule. 37 
The operator V c (re) is the crystal field potential energy, where Rp and —Zpe (<0) are the 
position vector and the electric charge of the i'th surrounding ion, respectively. As to A of 
eq. (7), C + (C_) is the coefficient in the case of an electron number of less than half (more 
than half) in an incompletely filled 3d shell. 35 The mean (r -3 ) is defined by eq. (8), where 
l^i) is the orbital state of the electrons in the cr-spin shell of the ground state (see eq. (41)). 
The quantity ^ B is the Bohr magneton, g is the g value, fio is the permeability, and h is the 
Planck constant h divided by 2tt. Also, S is the magnitude of S, and HL is the total orbital 
angular momentum. 

2.2 Vibration displacement 

For the Hamiltonian in the equilibrium state of §2.1, we take the effect of the vibration 
displacement into consideration. Regarding the vibration, we consider a simple model, in which 
only a single magnetic ion vibrates, while the surrounding ions are rigid. This model may be 
relevant to, for example, a single magnetic ion surrounded by several ions in the substrate. 
Here, the surrounding ions are assumed to behave as an oscillator with the magnitude of 
vibration displacement, which is much smaller than that of the magnetic ion. This assumption 
is justified under the following conditions: (i) The coupling within the surrounding ions is 
stronger than that between the magnetic ion and the surrounding ions, (ii) The total mass of 
the surrounding ions is sufficiently larger than the mass of the magnetic ion. 38 ' 39 For example, 
in the case of the Fe ion on the CuN surface, 10 the Fe ion is surrounded by two N ions and 
five Cu ions. The total mass of the surrounding ions is more than six times as large as the 
mass of the Fe ion. (iii) A force constant between the surrounding ions and the other ions 
adjacent to them in the substrate is much larger than that between the surrounding ions and 
the magnetic ion. 38 

We now focus on the effective nucleus of the magnetic ion. Here, the nucleus consists of 
a nucleus and core electrons, where the core electrons correspond to the electrons other than 
those in the u-spin shell. The effective nucleus is first assumed to have a slight vibration 



5/55 



J. Phys. Soc. Jpn. Full Paper 

displacement, Ar n (see Fig. 1 and Appendix ), i.e., 

Ar n = (Ax n ,Ay Q ,Az Q ). (10) 

The displacement of the £th electron in the u-spin shell, A77, is then given by 

Ar £ = (1 - ri)Ar n , (11) 

where — r?Ar n represents the difference in displacement between the nucleus and the electron. 
This —i]Ar n appears in the presence of the crystalline electric field due to the surrounding ions 
(see Appendix ). Here, the dimensionless quantity characterizing the difference r] (0 < r] < 1) 
has been assumed to be independent of £ and Ar n for simplicity. For example, this r\ may be 
written by eq. (A-8). Regarding eq. (11), we have considered the following situation: First, 
the effective nucleus of the magnetic ion has the displacement of Ar n . Second, each electron 
moves in a circular orbit around "the nucleus fixed at Ar n " at a high speed, and simultaneously 
experiences the crystalline electric field due to the surrounding ions. This electric field shifts 
the center position of the circular orbit from Ar n to (1 — r/)Ar n (see Appendix ). 

Thus, when the nucleus has the displacement of Ar n , the position of the £th electron 
measured from the nucleus, r' f , is represented by 

r' e = r e + Ar e - Ar n = r e - r/Ar n . (12) 

From the time differential of eq. (12), we obtain the following momentum: 

Pe = Pf ~ V^Pc (13) 

with 

Pe = mr' e , (14) 

Ap e = (Ap e , x , Ap e>y , Ap E)Z ) = mAr n . (15) 



2.3 Unperturbed and perturbed Hamiltonians 



We next derive the unperturbed and perturbed Hamiltonians from a Hamiltonian, % dlsp , 
which contains the effect of the above-mentioned vibration displacement. This % dls P can be 
obtained on the basis of % of eq. (1). Namely, when r' £ of eq. (12) and p' £ of eq. (13) are taken 
into account, H is rewritten as 

n disp = nf sp + ^ isp , 



N 



7^ + V en (r' e ) + V c (r' e + Ar n ) 



(16) 
(17) 

(18) 
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(19) 



r' e = y/(x e - r/Ax n ) 2 + (y e - r]Ay n ) 2 + {z e - 7]Az n ) 2 , (20) 

r' e + Ar n = r e + (1 - r?)Ar n . (21) 

For eq. (16), we expand V en (^) =^en(^-??Ar n ) and V c (r' g +Ar n ) =V c (re+(l-r})Ar n ) in 
power series in r/A^ n and (1 — r/)A^ n , respectively, with £=cc, y, z. In addition, (r' 3 ) of V^ sp 
is expanded in power series in rjA^ n . Since |?yA£ n |, |(1 — 77)A£ n |, and |r7Ap e ^| are considered 
to be sufficiently small, we neglect terms higher than the first order of r/A£ n , (1 — r/)A^ n , and 

For such % dlsp , we adopt Van Vleck's approach; 23 that is, the crystal field potential en- 
ergy in the equilibrium state is contained in the unperturbed Hamiltonian, while the spin-orbit 
interaction and the modulation of crystal field potential energy due to the vibration displace- 
ment are treated as perturbations. We therefore present the following effective Hamiltonian 
consisting of the unperturbed term T-Lq and the perturbed term V p : 

n cS = H + V P , (22) 

V p = V; + Vl n + V^ + V so , (23) 

^SO = ^ S o0 + ^sol + V so2 , (24) 



with 



^o = E 



N r 2 

Pe 



e=i 



^ + V en (re) + V c (r e ) 
2m 



vl 



m 

N 



Vi n = -„g (Ax n A + A, n A + Az n A) VWr<), 

VL = (l-^)g + A ^ + ^), 

V S o0 = AX • 5, 

^sol = »/C± ( -^Ax n + -^Ay D + -^Az n )L-S, 
V S0 2 = -??A- (r t x Ap c + Ar n x p t ) ■ S, 



where 



N 



rt = {x t ,y t ,z t ) = J2 r ^ 



(25) 
(26) 
(27) 

(28) 

(29) 
(30) 

(31) 
(32) 
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N 

Pt = (Pt,x,Pt,y,Pt,z) = J^P/, (33) 

e=i 

(l A ^) = £ (*i |3^ 5 A^ n | / f> , (34) 
x 1 e=i I e=i 

for £=x, y, 2. Here, eq. (25) is the zeroth-order term of r]Ap e ^, r]A^ n , and (1— rf)A£ n in eq. (17). 

Equation (26) is the first-order term of r]Ap e in YleLi P'l '/^m, while eq. (27) is that of f]A^ n in 

Y$Li v en(r' e )- Equation (28) is the first-order term of (1 - r/)A£ n in Y^=i v c( r i + (1 ~ »?)Ar„). 

In addition, eq. (29) is the zeroth-order term of 7/A£ n and i]Ap c ^ in eq. (18), while eqs. (30) 

and (31) are the first-order terms of them in eq. (18). Here, eq. (30) (eq. (31)) corresponds 

to a term with a change of A due to the vibration (a term with a change of L due to the 

vibration). Note that Ax n , Ay n , and Az n in ( ) of eq. (30) can be put outside of ( ), which is 

related to electron systems. Also, the suffix t in eqs. (32) and (33) represents the total sum. 



2.4 Orbital and spin states 



For the calculation of the perturbation energies, we here propose a specific orbital and 
spin state for %q of eq. (25). Simultaneously, we introduce a 3d orbital state that can be 
utilized in the principal axis transformation of §3.3. 

On the basis of eq. (25), we first have the 3d orbital and spin state with the largest S, i.e., 

\*j)\S z ), (35) 

where \S Z ) is a spin state with a (25 + l)-fold degeneracy of S z =— S - S, and |<3?j) is an orbital 
state of the electrons in the partially filled cr-spin shell. When the number of the electrons is 
N, the wave function for <&j, is given by the following Slater determinant: 

da(l) d a (2) ■■■ d a (N) 

d p (l) d p {2) ••• d p (N) 



(36) 



dy(l) d 1 (2) ■■■ dy(N) 
The suffix j of <l?j specifies the electron configuration, in which the occupied orbitals are 
represented by a, f3, ■ ■ ■ , and 7. The function d^(£) represents the 3d orbital of the <7-spin, 
where £=1, 2, • • • , N and (=a, /?,•••, and 7. We also have the following eigenvalue equation: 



Ho\* j ) = e j \$ j ), 



with 



= £d a + £dp + 

h 2 



-^Vi + V en {r £ ) + V c (r e ) 



d c (£) = s d( d c (£), 



(37) 

(38) 
(39) 
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where ej is the eigenenergy for and Ed ( is that for d^(£). In addition, |<E>i) is assumed to 
be a nondegenerate ground state with e\. Since "Ho is a real function, <E>i can be a real one. 
This \$>j) will be used to perform the principal axis transformation as described in §3.3. 

Second, we propose the orbital and spin state to obtain the perturbation energy (i.e., the 
spin-atomic vibration interaction and the spin- flip Hamiltonian) ; that is, 

\*j)\Sz), (40) 

where \S Z ) is the above-mentioned spin state, and \^>j) is an orbital state of the electrons in 
the <7-spin shell. The wave function for , is assumed to be 



Hsn. 



^a(l) 

M 1 ) 



^ a (2) 

M 2 ) 



^ 7 (1) V 7 (2) 



(41) 



with 



5 

E 

m=l(#C) 



^ = c ( kw+ E 4%w+E4 c >^) 



n=l 
-1/2 



1+ E 



,(0|2 



(C)|2 



(42) 



(43) 



where iV is the number of the electrons. The suffix j of ^fj specifies the electron configuration, 
in which the occupied orbitals are represented by a, f3, ■ ■ ■ , and 7. The function d^{£) is the 
dominant 3d orbital of the cr-spin, while d m {£) and p^(£) are the other 3d orbital of the cr-spin 
and the 4p orbital of the cr-spin in the same magnetic ion, respectively. Owing to the d-d 
and d-p hybridizations in the magnetic ion, d m (£) and p^(£) are included in eq. (42). The 
hybridization originate from, for example, the mixing of atomic orbitals via the surrounding 
ions. The bars of d m (£) and p^(£) are a mark to distinguish the components hybridized to 
d^{£). The quantities and Cp} are the coefficients for d m (£) and p^(£), respectively, and 



is the normalization factor, where |c^| 2 <Cl and \c y pj\' 2 <Cl are assumed. Note that d m (£) 
and p^(£) play important roles in the presence of the spin-atomic vibration interaction, as will 
be described later. 

We will also use the energy of in the perturbation calculation. For simplicity, the 
energy is assumed to be 



(C)|2 



Ej = e^iwoitfj), 



(44) 



where it is noted that Uq of eq. (25) (or eq. (2)) is the Hamiltonian in the subspace of the 3d 

E j = e j + Ae j , (45) 



orbitals. We here write Ej as 
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where ej is eq. (38), while Aej is given by Ej — ej. This Aej consists of terms higher than the 
first order of cjp and of] , where \Aei — Aej\ <C \&i — &j\ for % 7^ j. 

In this paper, eq. (44) is simultaneously regarded as the matrix element of the effective 
Hamiltonian of "Ho- Namely, we neglect the off-diagonal matrix element (^fji\7lo\^fj) for j 7^ j', 
which appears due to the hybridization, because the magnitudes of the off-diagonal elements 
are much smaller than \Ej — Ey\. The effective Hamiltonian will be used in eqs. (59) and (60). 

Furthermore, |^i) is assumed to be a nondegenerate ground state with E\. The function 
^1 is regarded as a real function because Hq is the real Hamiltonian. 40 

3. Spin- Atomic Vibration Interaction and Spin-Flip Hamiltonian 

We obtain the first- and second-order perturbation energies of the model in §2. By means 
of the principal axis transformation, we finally derive the spin-atomic vibration interaction 
and the spin-flip Hamiltonian as well as the conventional anisotropy spin Hamiltonian. 



3.1 First- order perturbation energy 



Using Vp of eq. (23) and l^)^) of eq. (40), we obtain the first-order perturbation energy 
i.e., 

V^ = (S z '\(^ 1 \V p \^ 1 )\S z ). (46) 

We here have V^j | x I / i)= (\Pi|V^|\I/i)=0, because L is quenched for the nondegenerate 
ground state whose wave function is the real function, i.e., (*i|X|*i)=0. 40-42 The first-order 
energy is then expressed as follows: 

vM^ + vg+vP+V™, (47) 

^S = (^'i(*ii^o2i*i)i^)= Yl u ^ ( 48 ) 



H=x,y,z 



= (S z '\^ l \V^ 1 )\S z 

N 



= -v E K 



£=x,y,z 







E^n(r,) 



#1 



2Mut 



V c {1) = (S Z '\^ 1 \V^ 1 )\S Z ) 

= (1-1/) £ (*, 



£,=x,y,z 







#1 



IMujj. 



( a ? + 4) 



Sz ,Sz ' 



V™ = <5,'|(*i|^|*i)|S,> 



E (*ibul*i) 1<; — 



= 0, 



(49) 



(50) 



(51) 



10/55 



J. Phys. Soc. Jpn. Full Paper 

with 

= -r/A- [(*i \yt\ *i) Ap c , z - (*i \z t \ *i) Ap 6)W ] 

= Ag] (-a, + 4) - Ag (-a, + 4) , (52) 
n y = -r/A^ [(*i \zt\ *i) Ap^ - (^i |x t | *x) Ap e , z ] 

= AW {-a x + 4) - Ag) {-a z + a\) , (53) 
IT = -t/A^ [(*i |x t | *i) Ap C)V - (*i *x) Ap e ,*] 

= AS + 4) " K vl (" a - + a *) > ( 54 ) 

where 



A U = _ ^ ft*M V 2""^ (55) 
In the last result of eq. (48), (S' z \ and \ S Z ) are omitted to use the operator representation. The 
vibration displacement of the effective nucleus, Ar n of eq. (10), has been replaced here with 

Ar n = {Ax n , Ay n , Az n ) 

= (^S: ( ax + 4 ) ' \f^ y ( ay + a ' vSi; ( az + a «)) • (56) 

The operator (a|) denotes the annihilation operator (creation operator) of the atomic 
vibration in the £ direction, 30 where the atomic vibration represents the vibration of the 
magnetic ion. The quantity M is the mass of the magnetic ion, and uj^ is the angular frequency 
in the £ direction of the magnetic ion. On the basis of eq. (15), Ap e has been written as 

Ap e = mAr n = m-^, (57) 

where Ap D is the momentum of the effective nucleus. This Ap n also has been replaced with 



= ^ (-ax + 4) , i^p {-a y + 4) , i^^i(-a 2 + 4)j . (58) 

Note here that Ar n (i.e., the displacement of the effective nucleus) may correspond to the 
position of the gravity point of the magnetic ion, because the mass of the nucleus is consider- 
ably larger than that of the electron. In addition, Ar n (Ap n ) has been assumed to be equal 
to the vibration displacement (momentum) of a mass point acting as a harmonic oscillator, 
in which a restoring force is linear in the displacement (see eq. (A-l) and Appendix ). The 
mass of the mass point is M, and its angular frequency in the £ direction is u^. 

In the derivation of eqs. (51) - (54), we have used (^ r i|pt i/1 |*i)=0 with fi=x, y, z. Here, 



11/55 



J. Phys. Soc. Jpn. 

ji\pt,ij\® j) is represented by 

TYl 

(*i'W*i) = ^(*j'lk,Wo]|*j) 
m(Ej - Eji) 



Full Paper 



-^Afit^. 



3h 



(59) 



where [fit, Ho] =J2eLi MiPt d Q m ) =1 ^Pt,^/ m and eq. (44) have been utilized. We also men- 
tion that IIu of eqs. (52) - (54) (or eq. (48)) is the ft component of (\I>i| — t?A^ (r t x Ap c ) 
where (\I/i|pt )jU |\l/i)=0 has been used in of eq. (48) with eq. (31). 

From now on, we will make no mention of eq. (49) or (50), which either contains only 
and but no spin operators. 



3.2 Second- order perturbation energy 

Using Vp of eq. (23), \^j)\S z ) of eq. (40), and Ej of eq. (44), we obtain the second-order 
perturbation energy , i.e., 

F (2) = y y ^l^ilUpl^OI^^/K^IUpl^i)!^) ^ (6Q) 

^ — ' ^ — ' Ei — Ej 

When terms higher than the first order of r/A£ n and rfAp e £ with £=x, y, z are neglected, 
is expressed as follows: 

K - v so0,so0 + v so0,sol + v so0,so2 + V so0,k' y 0L > 

^o0,so0 = /(^soO^soo) = E A^S,, (62) 

/j,,v=x,y,z 

^s ( od,sol = /(^o0, ^sol) + /(Usol, V so0 ) 
^so0,so2 = /(^soO, ^so2) + /(^so2, ^soo) 

= E (64) 

H,v=x,y,z 

V^ ok = f(V so0 X) + f(V^V so0 )= £ fc^, (65) 



H=x,y,z 



where 



F) = V V (^K^il^l^OI^X^/K^-iy^i)!^) . (66) 

^ — ' ^ — ' Ei — Ej 

S z " 3 

In the last results of eqs. (62) - (65), (S' z \ and \S Z ) are omitted to use the operator represen- 
tation. The quantities A^, A^ jI/; g, and fe M are given by 

A^ = A 2 £ M)1 „ (67) 
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/<■« Z//A " (r-3) \/2M^^' 
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r X)X = (-A+ ZiJ/ + A x . Z! y)(a y + a\) + (A+ ;j/)2 - A w )(a 2 + 4), (69) 

^x,y = ^x;z,x a x ~ ^-x;z,x a x + ^-y;z,y a y ~ ^y;z,y a \ + (~^x;x,z ~ ^y;y,z) a z + (^y-y,z + A :r;:r,z)4i 

(70) 

r^.z = — A x . yx a x + A^.^^aJ, + (A x . x ^ y + A z . z ^ y )a y + (— A^ ZiJ/ — A x . Xjy )a y — A z . y ^ z a z + A^. y ^ z a\, 

(71) 

^y,x = ~A X ; Z ^ x a x + A^.^^aJ, — Ay. zy a y + A y . zy ay + (A y .^ y ^ z + A,,.. 3,^)0,2 + (— A x . x ^ z — A y . yz )a z , 

(72) 

r s/,j/ = (K,z,x - K;z,x)( a x + 4) + (" A ^ + A ^,J(a 2 + 4). ( 73 ) 
= ( — ^-y,y,x ~ A z - tZiX )a x + (A^^ + A^.^^J, + Ay. xy a y — Ay. xy ay + A^.^^a^ — A^.^^a^, 



(74) 



F^rr — A^.^Oj; A+^4 + ( A z^z,j/ ^x;x,y) a y + + A z;z,j/)4 + ^z;y,z a z A^, Jl 



t 

y^zi 



(75) 



— (A^^ + Ay. yx )a x + ( Ay. y ^ x A ., . ., ^ .j, ) 4 A y . j. > y a y + A^^a^ A^t^a^ + A 2 ,. a , ; ,a 



t 



r z ,z — ( A^^, + A z . J/j;E )(a a ; + 4) + (At; X ,y A z . xy )(a y + 4)) 

i=x,y,z 

respectively, where 



£i - Ej 



(*i'l^l*i) = i 



m Mhuj£ Km 
MV 2 y 2Mw e ¥ 



(*j'H*i)> 



p = i_ ,, 

M V 2 



2Mw ? 



(76) 
(77) 
(78) 

(79) 
(80) 
(81) 

(82) 
(83) 
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In the derivation of eq. (80), we have used eq. (59), (^j'\L\^j) = — (^j\L\^ji), 42 and 
(^j'\Ct\^j) = (^j\£,t\^j')- In addition, eq. (82) has been obtained by using eqs. (83) and 
(59). Here, we mention the relations of = (^j\(F^\^ jf ) = and 

ji\F^\^ j) = (^j\F^\^fji), where {^j,\u t \^j) = (ty j\v t \%> . We also show the relations of 
= £»,», 43 ^,u=A u ,„ A,,,. .\,,„ :£ . A^ ;; , t A ( ; i:; , =(AJ.^)*, 44 and (T^=T v>lt , where 

Note that the two terms consisting of V^ o and V^ n + V£ in eq. (60) cancel each other, i.e., 

/(^soO, C + K) + /(Kn + K', ^so0) = 0, (85) 
where f(X,Y) is given by eq. (66). Here, we use (^j\L\^i)= — (\E'i|.L|\E'. / ) (see ref. 42) and 
(Wn + K'l*i)= (*ilKn + K'l^)- In particular, /(VW), t£) + f(V^V so0 )=0 corresponds 
to the so-called Van Vleck cancellation. 23,45 



3.3 Principal axis transformation 



To obtain expressions of °f e Q- (48) and of eq. (61) in the principal axis coor- 
dinate system, we first introduce the rotational transformation matrix for the principal axis 
transformation. This transformation is based on the conventional theory for the anisotropy 
spin Hamiltonian. 40 ' 46 Here, x, y, and z are regarded as the initial axes. 



In accordance with refs. 40 and 46, we focus on ^ 

,(2) 



Cn. vS,,S v included in V r 



(2) 



H,v=x,y,z '~ii,v>->^v 



so0,so0 



of eq. (62) and V^q so1 of eq. (63). The expression in the matrix representation is written as 





( c 

*^x,x 


r 


r 

*~-x,z 


\ 


( 


Sx 


\ 


] £n,vS = (s x , Sy, s z ) 


£y,x 


£y,y 


r 






s 

y 




H,u=x,y,z 


y C-z,x 


£z,y 




) 


V 


Sz 


J 



= SCS T , 



(86) 



with 





( c 


r 

x,y 


r 

'-'X^Z 


\ 


c = 


r 

y,x 


r 

^y,y 


r 

>~y,z 






y £z,x 


£z,y 


C-z,z 


J 



£(0) 



/ AO) AO) 

'-X'X 



V 



AO) 
AO) 



x,y 

(0) 
y,y 

(0) 
z,y 



£ (0) \ 

*"X,Z 

c 



,(0) 

^y,z 

,(0) 



(87) 
(88) 
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( £ xx £'„,, £,, \ 



£(0) 



J x,z 



V,x ~y,y ~y,z 

\ £ £ £ I 



ei 



r< - r - r (°) 



(89) 



(90) 



(91) 



where S 1 is the transposed vector of S=(S X , S y ,S z ). Here, is described by <E>j consisting 
of only d^(£), whereas £ is related to d m {l), p^(£), and d^{£). We also have £^°l=.ci ^ and 

/*' =/*' 43 

The principal axis of the crystal, named X, Y, and Z, are considered to be the coordinate 
system that diagonalizes £(°). 40 > 46 We then introduce the matrices U and U T that diagonalize 
C^°\ Here, U is the rotational transformation matrix, while the transposed matrix of U, U T , 
satisfies U T =U~ 1 . We thus have 

\ 



UC^U T = 



( £ {0) 







,(0) 

-Y,Y 









( c 



UC'U T 



x,x 



£ 



XY 



C 



Y,X 



£ (o) , 

£ \ 

C 



Y,Z 



\ C 'z,x 



£ 



Z,Y 



£ 



(92) 



(93) 



z,z ) 



where 



U 



( COS X ,x 

cos e Y , x 



u 1 



AO) . 
4,J 



COS X ,y COS X ,z ^ 

COS 0y,y COS 6y,z 

\ COS 9 Z ,x COS 9 Z ,y cos 6 Z ,z J 

( cos6 x ,x cos 0y,x cos6 z , x ^ 

COS 9 X ,y cos 0y ;J/ COS 6»Z,y 
y COS#x,z COS 9y,z cos 9z,z / 

cos 9 Itlx cos 1, „£$},, 



E 

H,v=x,y,z 



E 

H,v=x,y,z 



cos 0/ ;jU cos 9 JtU £ fll/ , 



(94) 

(95) 

(96) 
(97) 



with /, J=X, Y, Z. The angle 0/ ;(U is the relative angle between the / and \i axes, with I=X, 
Y, Z and /x=x, y, z. By utilizing Z^ ^/^ 43 £ I3 = Y.^ v=x ^ z cos9j, v cos9 I ^£ v ^ =£ }1 is 
obtained. 
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^-n,uS^S u — (S x , Sy, S Z )U T U 



ix,v=x,y,z 



( c 

£y,x 
\ ^z,x 



c c \ 



-y,y 



-z,y 



y,z 



( s x \ 



U T U 



= (Sx, Sy, Sz) 



( £ {0) 





Y,Y 








£ (o) . 



J 

A ( Sx \ 

Sy 
\ S z I 



Sy 





( 


^x,x 


C 

^X,Y 


^X,Z 


\ 


( 


Sx 


\ 


+(Sx, Sy, Sz) 




C 


c 

Y,Y 


c 






S Y 






V 


c 

^Z,X 


^Z,Y 


c 

^z^z 


) 


V 


S z 


) 



(98) 



= E 4°K 2 + E ^.jSjSj, 

I=X,Y,Z I,J=X,Y,Z 

where the relation between (Sx, Sy, Sz) in the principal axis (X, Y, Z) and (S x , S y , S z ) in 
the initial axes (x, y, z) has been defined by 



(99) 



3.4 Expressions of spin- atomic vibration interaction and spin- flip Hamiltonian 

Using the above-mentioned U and U T , we can obtain the expressions of of eq. (48) and 

Vi2) (= y soO,soO + ^s ( oO,sol + F s ( oO,so2 + V So^ ° f eq ' (61) in the P rinci P al axis coordinate system. 
Each term is named any one of the spin-atomic vibration interaction, spin-flip Hamiltonian, 
or anisotropy spin Hamiltonian. In addition, we consider the feature of the respective terms. 
Their coefficients and operators are listed in Table I. 



( Sx\ 




( 


Sx 


\ 


Sy 


= u 




Sy 




\ S z J 




\ 


s z 


J 



3 -4 A v soO,soO,l and u so0,so0,2 



Using eq. (98) , we first rewrite U go Q go0 of eq. (62) as 



u 



(2) 



(2) 

y soO,soO,l 



.(2) 

-^so0,so0,2' 



soO,soO 

^oO,soO,i = D S 2 z + E(S\ - S Y ) + FS(S + 1), 
-£U 2 = E ^ZSrSj, 



with 



A 



,(0) 



-z,z 



(o) , AO) 

X,X Y,Y 



D = 

E = X\C^ x -C^ Y )/2, 
F = A 2 (£^ + 4°) y )/2, 



)/2 



(100) 
(101) 
(102) 

(103) 
(104) 
(105) 
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Afj = A 2 £' /iJ; (106) 

where we have Af^=Aj^ using C' IJ =C I JI (see §3.3). The Hamiltonian ^goOgoOi is the so- 
called anisotropy spin Hamiltonian, 40 where D is an anisotropy constant. The Hamiltonian 

(2) 

u so0so0 2' wn i°h contains the spin-flip operators such as SxSz and SySz but no atomic 
vibration operators, is named the spin-flip Hamiltonian. 

From now on, FS(S + 1) in eq. (101) will be neglected because of a constant term. 



34-2 1A and 



( 2 ) 

Using eq. (98), we express V^ q q so1 of eq. (63) as 

T/ (2) _ (2) (2) 
*so0,sol u so0,sol,l u so0,sol,2> 



with 



v {2) 
so0,sol, 



1= [DS 2 Z + E(S X -S Y ) + FS(S + 1)] £ G^ + a\), 



€=x,y,z 



V 



,(2) 

so0,sol,2 



E E Xi,j£SiSj{at + al), 



£=x,y,zI,J=X,Y,X 



(r- 3 ) V 2Mu{ 



A'/../ :£ = 2r?A 2 



h 



(107) 
(108) 

(109) 



(110) 
(111) 



where we have A' 7 j.^=Aj^ using £J 1J =£J J1 (see §3.3). The interactions ^0 sol l an< ^ 

(2) ' 

Vg Q Q sol 2 , which contain the spin and vibration operators, are named the spin-atomic vibration 
interactions. 



3-4-3 VX 2 



Using U and U 1 , we rewrite V^q so2 of eq. (64) as 



(2) 



so0,so2 



E r M^^^ 



H,i/=x,y,z 





( 


r 

1 x,x 


r 


r 






( 




\ 


(S x ,S y ,S z )U T U 




r 

1 y,x 


L y,y 


r 




U T U 




Sy 






{ 


r 

z,x 


1 z,y 








\ 


s z 


J 



^ ^X,X Fx,Y ^x,Z ^ 

(S X ,S Y ,Sz) T Y ,X ?Y,Y T Y ,Z 

\ ^Z,X ^Z,Y ^Z,Z j 

E r i,J s i s Ji 



( Sx \ 

S Y 
V 5z j 



(112) 
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(113) 



with being eqs. (69) - (77). We here mention a relation of (Tj y 

52w= x ,y,z cos Oj," cos ^/,M r ^= r J,/ usin S ( r M,^) t=r ^M ( see § 3 - 2 )- Tne interaction I '^ )0 
which contains the spin and vibration operators, is the spin-atomic vibration interaction. 



) f = 

,(2) 



344 v£l k 



(2) 

We write k of eq. (65) by 



V 



(2) _ V- u a 
soO,k ~ 

H=x,y,z 



{k x , ky,k z )U U 



Sy 



with 



fc/ = ^ fc M cos6»/ )At , 



(114) 



(115) 



jj,=x,y,z 



(2) 

where A;^ is eq. (78). The interaction V^q k , which contains the spin and vibration operators, 
is the spin-atomic vibration interaction. 

8.4.5 V« 

We express °f e Q- (48) by 



H=x,y,z 



( s x \ 



= (U x ,U y ,U z )U T U 



J y 

\ s, ) 

= U X S X + TlyS Y + U Z S Z 



with 



IT/ = ^ n M cos6»/ iM , 

H=x,y,z 



(116) 



(117) 



where IT X , II y , and II 2 are eqs. (52) - (54), respectively. The interaction V^, which contains 
the spin and vibration operators, is the spin-atomic vibration interaction. 
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3.4-6 Considerations 

On the basis of §3.4.1 - §3.4.5, we finally express the perturbation energy, + , as 

v - v w +v w + v w W 2) W 2) +v {2) +v {2) 

v — v so2 + u soO,soO,l + u so0,so0,2 + u soO,sol,l + u so0,sol,2 + 1 so0,so2 + 'soO.k 

= Va + V SA + Vsf, (118) 

where Va., Vsa 5 and Vsf are the anisotropy spin Hamiltonian, spin-atomic vibration interac- 
tion, and spin-flip Hamiltonian, respectively. These terms are written by the following simpli- 
fied expressions: 

Va = v { XoO,l = DSl + E(S 2 X - S 2 Y ), (119) 

V SA - V so2 "I" u so O,sol,l + °so0,sol,2 + ^so0,so2 + "soO,k 



= E E c /,5^(-«« + 4) 

i=X,Y,Z£=:r,y,z 



/,j=x,y,zf=x,j/,« 

+ [£S 2 +£(Si-S 2 ) + FS(S + l)] ^ G^ + a\), (120) 

^-^300,2= E A Z S ^ ( 121 ) 
I,J=X,Y,Z 

where c/^ is the coefficient of the Sj(—a^ + a^) term, and c\j t j£ (c2,i t j^) is that of the SjSja^ 
term (that of the SiSjat term). Here, c/^ consists of of eq. (55) and A*^ of eq. (80), 
while Cnj,j£ (n=l, 2) is composed of A'j j.^ of eq. (Ill) and A^.^ of eq. (81). In other words, 
V s ( 2 and V s ( o 2 q k are terms with Si(-a^ + a^), while v^ Q sol 2 and V s ( q 2 q so2 are those with SiSja^ 
and SjSja^. 

We consider the features of V g % V®^, «S) )S ol,i' v soO, S ol,2> y S oO,so2 of y SA, and ^00,2 

(1) (2) (2) (2) 

of V SF . We first find that V g ^ and V^ k are proportional to A, while </ 0)Soljl , </ ,sol,2> 
^so0so2> anc ^ u so0so0 2 are P ro P rti on al to A 2 (see Table I). Second, the presence or absence 
of the respective terms is shown by examining the respective coefficients for each set of rj, c p , 
Cd, where and c p are defined as c^=c^ and c p =c p ^ , respectively (see Table I). Here, lo^ is 
put to be oj^=2irf^ for £=x, y, z, where /e is the vibration frequency in the £ direction. In 
addition, e\ — ej=E\ — Ej is roughly set by taking into account |Aei — Aej|<C|ei — ej\ for 
j ^ 1 (see eq. (45)). The details are described as follows: 



Spin-atomic vibration interactions Vsa 
Sj(—a^ + a|) term 

- V™ with A« oc r?^ (see eq. (55)) 

The interaction V^ appears for rj^O and c p /0, because each term contains 
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ri(^fji\/j,t\^j), with fj,=x, y, z. Here, (tyj'\(it\^j) is the same formula as a matrix el- 
ement of the so-called electric dipole transition; it is nonzero for c p ^0. In addition, 
is proportional to y/fe, reflecting that contains Ap n ^ of eq. (58), with £=x, y, z. 

y soO,k With A M,« K ^Vh ( See e( l- ( 8 °)) 
(2) 

The interaction K, g k , in which each term contains 77 ( ^ ^ | I ^ 1 ) with £=x, y, z, is 
present for rj^O and c p /0. In addition, V^j k is proportional to because V^j k (or 
VQ contains Ap n ^. 

SiSja^ and SiSja^ terms 

- y so0,so2 with K %v£ « »7 ( see e q- ( 81 )) 

(2) 

The interaction V^ q q so2 appears for and c p ^0, because each term has 7](^j/\u t \^>i), 
with i/=x, y, z. In addition, V^j go2 contains both terms with yj~f^ and those with 
l /\fh, reflecting that V^ so2 (or V so 2) has Ap n ^ and A£ n of eq. (56), with £=x, y, z. 

«S),sol,i with G ? K ^r' 5 )/^ (see eq. (110)) 

The interaction gol 1 exists for 7/7^0 and c p j^0, because it is proportional to r/(3£r~ 5 ), 
with £=x, y, z. Here, (3£r -5 ) is nonzero for c p ^0. It is noted that the operator is 
proportional to £ as with the matrix element of the electric dipole transition. In addition, 
v so0,sol,i is proportional to l/y/fe because Vgo0,sol,i ( or ^sol) contains A£ n . 

- «S),sol,2 with A 'l,J;Z K (s^ eq. (Ill)) 

(2) 

The interaction v so q gol 2 is present for rjf^O, c p ^0, and c^O, because it is proportional 
to 7]{?>^r~ 5 )C'j j with £=x, y, z. Here, C'j j is nonzero for c^O. In addition, so1 2 is 
proportional to 1/y^, reflecting that gol 2 (or Vgoi) has A£ n . 

• Spin-flip Hamiltonian Vsf 

- ^o0,so0,2 with A Z * £'l,J (see eq. (106)) 

The Hamiltonian 02 , which is proportional to C'j j, exists for q /0. Note that 

(2) ' ' 

v soQ soO 2 ^ S independent of the atomic vibration, and it does not contain rj and f%. 

4. Application to Fe Ion 

As an application, we consider the spin-atomic vibration interaction and the spin-flip 
Hamiltonian of the Fe ion (Fe 2+ ) in a crystal field of tetragonal symmetry, where Fe 2+ has 
six 3d electrons (i.e., 3d 6 ). We first investigate the / dependence of their coefficients. Second, 
their coefficients are evaluated by assigning appropriate values to A. 

4-1 Orbital and spin states 

As to the Fe ion (3d 6 ), the ground state consists of five up-spin electrons and one down- 
spin electron, according to Hund's first rule. The spin S is therefore considered to be 5=2. In 
addition, since the up-spin shell is filled, *j of eq. (41) is given by the orbital wave function 
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A, E 5 (E 3z2 _ r2 ) 

K E 4 (E x 2_ y2 ) 

A ° — E 2 =E 3 (E x =E yz ) 
E l (E xy ) 

Fig. 2. Energy level Ej of the down-spin shell of the 3d orbitals of the Fe 2+ (3d 6 ) in the crystal 
field of tetragonal symmetry. We use the notation of Ei=E xy , E 2 —E xz , E 3 =E yZ7 E A =E x 2_ y i 1 
and E 5 =E 3z 2_ r 2, where xy, xz, yz, x 2 — y 2 , and 3z 2 — r 2 correspond to d\, d 2 , cJ 3 , c£ 4 , and d$, 
respectively (see eqs. (125) - (129)). Here, the relation of E x < E 2 (= E 3 ) < E± < E 5 is considered, 
where the second excited states are doubly degenerate, i.e., E 2 =Ej,. In addition, A s , A TO , and A; 
represent E 2 — E\ , E4 — Ei, and E$ — E\ , respectively. 



of the down-spin electron, i.e., 



v ", <•, c ,['!,■ E 4^m+E c 2^h ( 122 ) 

m=\{+j) n=l J 

^=( 1 + E i^i 2 + Ei«i 2 • ( 123 ) 

m=l(^j) n=l I 



Equation (36) also becomes 

$j = dj. (124) 

The energy of is given by Ej = ej + Ae-,- (see eq. (45)), while that of is ej (see eq. (38)). 
Here, Ej is assumed to be Ej in Fig. 2 by taking into account the crystal field of tetragonal 
symmetry. The wave functions dj and pj are then written as 

di = xyf 3t2 (r), (125) 

d2 = xzf 3 , 2 (r), (126) 



d 3 = yzf 3:2 (r), (127) 
1 

2 1 



di = \{«? ~ y 2 )h,l(r), (128) 



and 



d 5 = ^=(3z 2 -r 2 )/ 3 , 2 (r), (129) 



Pi = s/4,i(r), (130) 

P2 = yf4,i(r), (131) 

P3 = 3/4,1 W, (132) 
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respectively. The function f n ,l( r ) with the principal quantum number n and the azimuthal 
quantum number I is given by 

M r ) = ~ J= m exp (-^—] , (133) 



/4,lW 



512\/57fa B {p 



20r / r 



80 + 



a B.p \ a B 



™ p I -^.1 • (134) 



OB.d = i^—, (135) 

^cff,d 

«B,p = TT^-, (136) 
^eff,p 

with a B =eo/i 2 /(vrme 2 ), where a B is the Bohr radius and -Z e ff,d (^eff,p) is the effective nuclear 
charge for the 3d electron (4p electron). 

4-2 Notation and parameter setup 

In this system, we first replace r t =(x t ,y t , z t ) and p t =(Pt,x,Pt,y,Pt,z) with (x,y,z) and 
(p x ,Py,Pz) of the one-electron system, respectively. The angular frequency uj^ is set to be 
=2-7r/, where / is a vibration frequency. For eqs. (122) and (123), we put c^—Cd 
for m=l - 5 with m / j, Cp}=c p for n=l - 3, and Cj=(l+4c^+3c^) -1 / 2 = C. In addition, we use 
the notation of ei=e xy , e2=e xz , e 3 =e yz , e^=e x 2_ y 2, e§=e 3z 2_ r 2, Ei=E xy , E2=E XZ , E 3 =E yz , 
E4=E x 2_ y 2, and E$=E 3z 2_ r 2. The energy difference is then expressed as E xz — E xy =E yz — 
E xy =A s (> 0), E x 2_ y 2 -E xy =A m (> 0), and E 3z 2_ r 2 -E xy =A t (> 0) (see Fig. 2). The energy 
difference e\ — e.j=E\ — Ej is roughly set by taking into account \Ae\ — Aej|<c|ei — ej\ for 
JVI- 

As a parameter set, we choose Cd=c p =0.15, 47 A s =0.45A m , A;=1.45A m , A m =0.5 eV, and 
/=0.01 THz - 100 THz, which correspond to typical atomic vibration frequencies. 48 The mass 
of Fe 2+ , M, is M=(26+30)x(1.67xl0- 27 )+24x(9.11xl0~ 31 ) kg, where the number of protons 
(neutrons) is 26 (30), the mass of the proton (neutron) is 1.67xl0 -27 kg (1.67xl0~ 27 kg), 
the number of electrons is 24, and the electron mass is 9.11 xlO -31 kg. The effective nuclear 
charge Z e g,d (Z c s,p) is evaluated to be Z e s^=6.2b (Z e ff )P =4.10) according to Slater's rules. 37 

Note that from now on we will use the coordinate system x,y, and z for some coefficients 
in accordance with the notation in §2. Actually, however, x=X, y=Y, and z=Z are realized 
in this system, as described in Appendix C.l. 

4-3 Frequency dependence of coefficients 

The coefficients of the spin-atomic vibration interaction are now written as A^l/X, A*; ,/A, 
G*£, F/D, A'j j.^/D, and A^.^/D, while that of the spin-flip Hamiltonian is written as A^ F j/D. 
These expressions are dimensionless quantities. Note that E of eq. (104) becomes zero owing to 
the tetragonal symmetry (see eq. (C-18)). The expressions of all the coefficients are described 
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/ (THz) 

Fig. 3. Vibration frequency / dependence of A^/(A?7i) for (j,, £= x, y, z. The parameters are set to 
be Cd~Cp— 0.15 and A m =0.5 eV. 

in Appendix . 

In Figs. 3 and 4, we show the / dependences of A^/(Aryi) and A^/(Aryi), respectively. 
These expressions are the coefficients of the Sj(—a^ + a|) term. The magnitudes of the coef- 
ficients increase with increasing /, because they are proportional to y/J. 

Figures 5-9 show the / dependences of A^ F j/D, G^/rj, A'j j.^/(Dr]), A^. u ^/(Drj), and 
A~ v ^/(Drj), respectively. These expressions are the coefficients of the SjSja^ and SjSja^ 
terms. The coefficient Af^/D, which contains no /, certainly takes a constant value. The 
coefficients G^/rj and |Aj j.^/(Drj)\, which are proportional to 1/y/J, decrease with increas- 
ing /. This behavior is in contrast to that of |A^/(Ar/i)| and |A^/(Ar/i)|. The coefficient 
\A^. u ^/(Drj)\ monotonically decreases with increasing / in a region of / < 20 THz, because 
the dominant term of re^ in A^.^/D (see eq. (C-80)) is ir](m/h 2 )y/h/(2Muj^)A u , where 
W£=27r/. In the vicinity of /=60 THz, \A~. U ^/D\ drops, reflecting that some k~^s approach 
zero. Here, k~ f =0 reduces to hujc=hf=A u . 

4-4 Evaluation of coefficients 

For the evaluation of the coefficients, we here consider a model in which the Fe ion (3d 6 ) 
in the crystal field of tetragonal symmetry vibrates in the x direction. The spin-orbit coupling 
constant A is assumed to be A=— 5, —10, —25, and —40 meV, where A of Fe 2+ has been 
previously evaluated to be —12.6 meV. 35 The anisotropy constant D can then be evaluated 
by substituting the above A's, Cd=c p =0.15, A m =0.5 eV, and A s =0.45A m into eq. (C-17). 
Namely, we have £>=-6.63x 10~ 2 meV for A=-5 meV, D=-2.65x 10" 1 meV for A=-10 
meV, D=— 1.66 meV for A=— 25 meV, and D=— 4.25 meV for A=— 40 meV. We also use 
77=0.05, which is evaluated in Appendix . 

Using the above A's and D's and the results in Figs. 3 - 9, we can evaluate the coefficients 
of Vsa and Vsf- In particular, we focus on large components in the respective coefficients, i.e., 
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Fig. 4. Vibration frequency / dependences of the coefficients. Upper panel: A^, l/ /(Xr)i) for i/=x, y, z. 
Middle panel: A^„/(A?7i) for v=x, y, z. Lower panel: ^(A^i) for v=x, y, z. Note that Ay X /(Xi]i) 
takes almost the same value as Ay y /(Xi]i) at each /. The parameters are set to be Cd=c p =0.15 
and A m =0.5 eV. 



A$, A* x , G x D, G x F, k' xz . x , k' YZ . xl and K% z , x in V SA , and kf z and kf z in V SF . In the 
upper and lower panels of Fig. 10, we show the / dependences of the coefficients with A=— 5 
and —10 meV, respectively. In addition, the upper and lower panels of Fig. 11 show those 
with A=— 25 and —40 meV, respectively. 

We first find that the magnitudes of the coefficients of Vsf are larger than those of the 
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Fig. 5. Vibration frequency / dependence of Af^/Z? for I, J= X, Y, Z, where Af^^Aj 1 ^ (see §3.4.1). 
The coefficients take the respective constant values. The parameters are set to be Cd=c p =0.15 and 
A m =0.5 eV. 
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Fig. 6. Vibration frequency / dependence of G^/r] for £=x, y, z. The parameters are set to be 
Cd=Cp=0.15 and A m =0.5 eV. 



Vsa- This relation mainly indicates that the coefficients of Vsa are proportional to r\ (i.e., the 
degree of the difference in vibration displacement), whereas Vsf does not contain 77. As found 
from eqs. (Ill) and (106), the major difference between A' £ j.^ and Af^ may be the presence 
or absence of r]. 

In Vsa, "\G x D\, \G x F\, \A' x z .J, and |Ay Z .J" tend to be larger than |A^.J and |A* X |, 
because of the relation of |V^ i| » | T^ Q 2 1 ; 1^1 (see eqs. (63) - (65)). In particular, the 
relation of |V^ i| » IV302I indicates that the rate of the change in A due to the vibra- 
tion, |AA/A|, is larger than that in L due to the vibration, \AL/L\, in a region with an 
atomic radius of r ~ 0.1 nm (see eqs. (30) and (31)). Here, AA and AL are given by AA= 
r]C±(dr~ 3 /dx)(-Ax n ) =r]C±(3xr~ 5 )Ax n and AL= -(r]/h)(r t x Ap c + Ar n x p t ), respec- 
tively, where Ap c =(Ap C:X , 0, 0) and Ar n =(Ax n , 0, 0) are set in this section. 
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Fig. 7. Vibration frequency / dependences of the coefficients. Upper panel: j. x /(Dr)) for /, J=X, 
Y, Z. Middle panel: A'j j. y /(Dr]) for J, J=X, F, Z. Lower panel: A' IJ:z /(Dr]) for /, J=X, F, Z. 
We have A' 7 j.^—A'jj.^ (see §3.4.2). The parameters are set to be Cd=c p =0.15 and A m =0.5 eV. 



5. Comparison between Vsa and the Conventional Spin-Phonon Interaction 

In this section, we compare Vsa with the conventional spin-phonon interactions derived 
by Mattuck and Strandberg 25 and by Hartmann-Boutron et al. 30 In particular, we focus on 
the / dependences and magnitudes of the coefficients in the case of the Fe ion. We also use 
7/=0.05, which is evaluated in Appendix . 
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Fig. 8. Vibration frequency / dependences of the coefficients. Upper panel: K^. u ^/{Dr]) for v, £=x, 
y, z. Middle panel: Ky. v ^/{Drj) for v, y, z. Lower panel: A* ^/(Dr)) for v, !;=x, y, z. We 

note that A+ ;a . ^/{Drj) takes almost the same value as ^/(Drj) at each /. The parameters are 
set to be Cd=c p =0.15 and A m =0.5 eV. 



5. 1 Modification of Vsa 

To compare Vsa with the conventional spin-phonon interaction (see §5.2), we first mod- 
ify Vsa to a novel spin-phonon interaction. We here consider the following one-dimensional 
lattice system: ions with the same mass M are placed at even intervals along the x axis, 
and vibrate in the x direction. In addition, a specific ion located at the origin possesses S, 
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Fig. 9. Vibration frequency / dependences of the coefficients. Upper panel: A~. v ^/(Drj) for v, £=x, 
y, z. Middle panel: A~ v ^/(Dr]) for u, £=x, y, z. Lower panel: A~. v ^/{Drf) for v, £=x, y, z. We 
note that A~ x ^/{Drj) takes almost the same value as A~. y ^/{Drj) at each /. The parameters are 
set to be Cd~c p ~0.15 and A m =0.5 eV. 



and also this ion is surrounded by several negative ions. The system may be just described 
by a simple model in which S located at the origin is coupled to the phonon for the lat- 
tice system. Thus, Vsa in this model is obtained by replacing Ax n in eq. (56) and Ap IhX in 
eq. (58) with (1/VN) Z q y/h/(2Mu q ) (a q + aj) and (1/ViV) Z q ^Mhw q /2 (-a q + aj) , re- 
spectively. Here, N is the number of unit cells, q is the wave vector, co q is the angular frequency 
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Fig. 10. Vibration frequency / dependences of A^, x , G X D, G X F, h! x z . x , Ay ZiI , and Af. z x in 
Vsa, and A|F z and A|- F z in Vsf- Upper panel: A=— 5 meV. Lower panel: A= — 10 meV. The cases of 
A=— 5 and —10 meV lead to D=— 6.63x 10 -2 and — 2.65xl0 _1 meV, respectively. The parameters 
are set to be 7y=0.05, Cd=c p =0.15, and A m =0.5 eV. 



of the phonon at q, and a q (a\) is the annihilation operator (creation operator) of the phonon 
with q. As a result, Vsa of eq. (120) is modified as follows: 

Vsa = Yl Sl ~7W Cl ' x 'i {~ a i + 4) 



'N 

i,j=x,y,z v g 
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Fig. 11. Vibration frequency / dependences of A^, A* G X D, G X F, A' x z . x , A.' Y z . x , and Af in 
VsA) and A|F Z and A|? z in Vsf- Upper panel: A=— 25 meV. Lower panel: A=— 40 meV. The cases 
of A=— 25 and —40 meV lead to D——1.66 and —4.25 meV, respectively. The parameters are set 
to be 77=0.05, Cd—c p — 0.15, and A m =0.5 eV. 



+ [DS 2 Z + FS(S + 1)] -_ £ G x , q {a q + a\). (137) 

ViV q 

The coefficients c I>x>q , c^i^g, c 2 ,i,j, x , q , and G x>q , respectively, correspond to c/ )X , c lt i t j !X , 
C2,i,j,x-, and G x in eq. (120), where uj x in these coefficients has been replaced with u q . 

We now replace oj q with 27r/ in order to investigate the / dependence of the coefficients. 
We then can regard cj jX , g , ci,/,j jXj g, c 2 ,i,j, x , q , and G x , g as c/ jX , ci,/,^, c 2 ,/,j, x , and G x in eq. 
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Fig. 12. Vibration frequency / dependence of the magnitude of coefficient of the spin-phonon inter- 
action, \c q \. Solid lines: The coefficients of Mattuck and Strandberg 25 estimated for |A|=5, 10, 25, 
and 40 meV. Dotted line: The coefficient of Hartmann-Boutron et al. 30 In this coefficient, A is 
unspecified. 



(120), respectively, where cu x =2irf. Such coefficients may be directly compared with those of 
the spin-phonon interaction of the below-mentioned eqs. (138), (139), and (140). In particular, 
we focus on large components in the respective coefficients in the case of the Fe ion (3d 6 ) in 
the crystal field of tetragonal symmetry, i.e., A^, A^ jX , G X D, G X F, A' x z . x , A.' YZ . X , and A.f. z ^ x 
(see Figs. 10 and 11). 

5.2 Spin-phonon interaction 

The spin-phonon interaction originally represents an interaction between a single spin 
located at the origin and the phonon of a lattice system. 25 ' 30 Here, a unit cell of the lattice 
consists of a magnetic ion (with a single spin) and the surrounding ions. In addition, all the 
ions vibrate. Note that exchange interactions between the spins are neglected. 

As a similar lattice system to that in §5.1, a one-dimensional lattice system is considered, 
in which each unit cell is aligned along the x axis. In addition, all the ions vibrate in the x 
direction. The dominant term in the spin-phonon interaction, 25 ' 30 Vsp, is then expressed as 

v SP = ^-^E( c ^ + <4)> (138) 

I,J=X,Y,Z ViV q 

where N is the number of unit cells, q is the wave vector, and c q is the coefficient. 
The coefficient of Vsp derived by Mattuck and Strandberg, 25 c q , is written as 

1 / hu q 6e 2 r 2 
° q ~ 4vr6oV 2M U *; 2 A2 J R3 A ' ^ 

with u) q =27rf q . Ag Here, A is the spin-orbit coupling constant of eq. (7), A is the energy splitting 
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under the crystal field, R is the distance between the magnetic ion and the surrounding ions, 
tq is the ionic radius, M u is the mass per unit cell, and v is the phonon velocity defined as 
v=uj q /\q\, where v is constant. This spin-phonon interaction is essentially the same as the ex- 
pression by Van Vleck, 23 which corresponds to the third-order perturbation energy. Here, the 
crystal field potential energy is the unperturbed term, while V^ o of eq. (29) and the modula- 
tion of crystal field potential energy due to the lattice vibration are the perturbed terms. 23 ' 29 
We mention that the second-order perturbation energy due to V^ D o and the modulation of 
crystal field potential energy vanishes owing to the Van Vleck cancellation. 23,45 

To estimate c q , we assume A=l eV, R=l nm, and ro=0.1 nm. We also use 
M u =[(26+30)+5x(29+34)+2x(7+7)]x(1.67xl0- 27 )+(26+5 x 29+2 x 7)x(9.11xl0~ 31 ) kg on 
the basis of results of the STM experiments for the Fe ion on the CuN surface. 10 The unit cell 
has been assumed here to consist of one Fe ion, five Cu ions, and two N ions. The numbers 
of protons (the number of neutrons) of Fe is 26 (30), that of protons (that of neutrons) of 
Cu is 29 (34), and that of protons (that of neutrons) of N is 7 (7). The numbers of electrons 
of Fe, Cu, and N neutral atoms are 26, 29, and 7, respectively. Since the numerical value of 
v is unknown for this system, we utilize f=1.45xl0 3 m/s, 50 which is the value evaluated for 



On the other hand, the coefficient of Vsp derived by Hartmann-Boutron et al. c q is given 



with u q =2ir f q . Here, q=uj q /v has been used, where v is constant. The quantity V q is re- 
garded as the coefficient with a unit of energy. 6 This interaction originally arises from the 
spin-orbit interaction and the local strain at the spin site due to the phonon. To estimate 
c q , we use | V^|=4.45 meV and the above-mentioned M u and v. This \V q \ is obtained from 
C pTe =3\Vq\ 2 /(27rh 4: pv 5 ), 1 with p=M u /a 3 , where the volume of the unit cell, a 3 , is set to be 
a 3 =l nm 3 . The quantity C prc corresponds to a prefactor of the expression of the relaxation 

rate, and then C pre =10 4 Hz/K 3 is chosen on the basis of experimental results for magnetic 

i ^1 
ions. ' 

5.3 Comparison 

We first make a comparison of the / dependence of the coefficients between Vsp and Vsa- 
Figure 12 shows the / dependence of \c q \ of eqs. (139) and (140), where / is defined by f=f q . 
We find that \c q \ is proportional to \ff . Such / dependence (i.e., y/f) is different from those 
of the coefficients of the SiSja q and SiSja\ terms of Vsa (i-e-, 1/y/J), but the same as that 
of the Si(—a q + a\) term of Vsa- 

Next, as seen from Figs. 10 - 12, \G X D\, \G X F\, \A.' X z . x \, and \A' YZ . X \ in Vsa are larger 
than \cq\ in the low / region. For example, when A=— 10, —25, and —40 meV, the former are 



Mn i2 . 



by 




(140) 
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larger than the latter in a region of / < 1 THz. Also, the region of /, in which the former 
become larger than the latter, expands with increasing A. 

As the main reason why the magnitudes of the coefficients of Vsa can be larger than \c q \ of 
Mattuck and Strandberg, we consider that M in the coefficients of Vsa (oc M -1 / 2 ) is smaller 
than M u in c q of Vsp (oc M„ -1 / 2 ). Originally, the decrease in M increases the magnitude 
of vibration displacement |A£ n | (£=x, y, z) (see eq. (56) and Appendix ). Furthermore, the 

(2) 

increase in |A£ n | enhances the magnitude of Vsa (i.e., v^Jq so1 2 ) ( e -S-i eo ± s - (63) and (30)). 



On the other hand, c q of Hartmann-Boutron et al. is rewritten as c q =i yjWojqV 3 C prc /(3a 3 ) , 
by substituting the above-mentioned | V q \ 2 =(2irh 4: /3)(M u /a 3 )v 5 C pre into eq. (140). Since the 
related quantities of V , & , £111(1 C^pre 

in c q are not contained in Vsa, it may be impossible to 
identify the cause of the larger magnitudes of the coefficients of Vsa than \c q \. 

We anticipate that at a low /, the 1/y/J terms in Vsa can be verified from the pressure 
effect of the spin relaxation time, and so on. 

6. Transition Probability per Unit Time 

In this section, we discuss transition probabilities per unit time due to Vsa and Vsf for a 
simple model of the Fe ion with bistable states. Our interests are in how the expressions of Vsa 
(e.g., SxSz(a + a^)) and Vsf (e.g., SxSz) influence the transition probabilities. Note also that 
the transition probabilities will be necessary for investigating spin relaxation phenomena. 19 

When Vj is considered to be a time-dependent perturbed Hamiltonian, we have the prob- 
ability per unit time of a transition from a state \u) to \v), W^ v : 

^ = |lW|u)|^ ^_ E ^ +(B/r)2 . (141) 

with Hq S \u)=E u \u), where T~Lq S is the unperturbed Hamiltonian for the atomic spin, and 
\u) is the eigenstate of T~Lq S . 52 Here, a correlation time for Vj, r, is introduced by assuming 
that the correlation function for Vj has an exponential decay. In addition, the case of E U =E V 
represents the conservation of energy for a transition from \u) to \v), while the case of E u / E v 
represents the nonconservation of energy for the transition. 
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Fig. 13. Energy levels of -\D\S% with 5=2. 
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We now propose Hq S and Vj of a simple model for the Fe ion in the crystal field of 
tetragonal symmetry. Here, this ion is considered to vibrate in the X direction. Note again 
that x=X, y=Y, and z=Z are realized in this system, as described in Appendix C.l. When 
the representative terms in Vsa of eq. (120) and Vsf of eq. (121) (see Figs. 3-8) are adopted, 
the Hamiltonian is written as follows: 

U$ s = DS 2 Z + frwa^a, (142) 

Vsai = A|5 z (-a + a t ), (143) 

V S A2 = [A' x ,z( S xSz + S Z S X ) + A' Y ,z(S Y S z + S Z S Y )} (a + a f ), (144) 

V SA3 = G [DS 2 Z + FS(S + 1)] (a + at), (145) 

V SF = Afz(S x S z + S Z S X ) + Af z (S Y S z + S Z S Y ), (146) 

where a (a T ) is the annihilation (creation) operator of the atomic vibration in the X direction, 
and uj (>0) is the angular frequency in this direction. The coefficients D, F, G, and A z are 
eqs. (CT7), (CT9), (C-62), and (C-29), respectively. In addition, A' x z and A' YZ correspond 
to eq. (C-76), while Af z and A| F Z to eq. (C-54). The state |u) is written as \S z ,n), where n 
is the vibrational quantum number. 

For this model, we choose S=2 (see §4) and D=—\D\ using examples from the experimental 
and theoretical studies of the Fe ion on the CuN surface. 10 As shown in Fig. 13, the model 
exhibits bistability between S z =— 2 and 2. The energy difference between the ground states 
and the first excited states is 3\D\, while that between the first excited states and the second 
excited state is \D\. 

By using eqs. (141) - (146), the respective W*_^„'s are obtained as follows: 

- >^H {r J+ W rr <147) 

Wl^, +1 = ||M|^(» + ^- (nwf % iT? , (148) 



2vr 



(A' ) 2 + {A'* 



W h stUs z ±i,n-i = T ~ 4 H*Sz ± l) 2 (S T S Z )(S ±S Z + l)n 

1 h/r 

X k[-\D\(±2S z + 1) -M 2 + (a/r) 2 ' 

(AV z ? + (AU) 2 



(149) 



W§^ Sz ±i,n + i = T l — 4 — J (2^z ± ir(S T S Z )(S ±S Z + I)(n + 1) 

x- ^ 5 , (150) 

Wli^-> = + + ^ (151, 
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Wl*Us z ,n + i = yG 2 [-\D\Sl + FS(S + 1)] 2 (n + 1)1^^1^, (152) 



<,n^ z± i, n = — , 1 ^(25 Z ± 1)*(S T 5 Z )(S ± S Z + 



2vr 



(A| F z) 2 + (A^) 2 



1 ft/r 



h 4 v ^ /v-t^a" > n D 2( ± 2S Z + I) 2 + (W 2 ' 

(153) 

where is an imaginary number, whereas the other coefficients are real number (see §4). The 
transition probability Wg^ 2 n ^ Sz±1 n _ 1 iWs z "n^S z ±i n+1 ) corresponds to a probability for a 
transition with the conservation of energy in the case of flu=—\D\ (±25^ + 1) (huj=\D\(±2Sz + 
1)). The other W^^'s are probabilities for transitions with the nonconservation of energy, 
irrespective of a; (> 0). 

The parameters are now set to be |D|=1.5 meV, 10 ?7=0.05 (see Appendix 
), A|=-(1.00xl(T 5 )|L>|i, A^ jZ =(7.02xlO -3 )|.D|, A^ z =-(1.94xlO -3 )|D|, G=9.40xl(T 3 , 
F=-1.25|£>|, A| F z =(7.47xlO- 1 )|D|, and Aff z =-(2.07xl0 _1 )|£>|. Here, A' xz , A' YZ , G, Af z , 
and Ay F z have been respectively evaluated from A' x z . x , A' YZ . X , G x , A xz , and A Y F Z at f=l 
THz in Figs. 5-7, while A k z has been roughly estimated from A k z x at f=l THz in the upper 
panel of Fig. 11. In addition, F has been evaluated from the above \D\ and eqs. (C-19) and 
(C-17) with A s =0.45A m . 

In Fig. 14, we show the Sz dependence of W^ v in the case of Hoj=3\D\. When t=1CT 5 s, 
i,n-i/ n > W^_ w /{n + 1), W?^ n+l /(n + 1), and W^^/n take the ex- 
tremely large values because of the transition with the conservation of energy. When r=10~ 10 
s, however, such probabilities get closer to W<j ^ n - > s z ±i n- ^ or eacn T i when a ratio in W^ v 
is defined as the largest W^ l _ yo divided by the smallest one, the ratio of r=10~ 5 s is larger 
than that of r=10 -10 s because of the sharpness of peaks of W^ v . Figure 15 shows the Sz 
dependence of W^_^ v in the case of hu=2\D\. All the W^_^ v 's are the probabilities for the 
transitions with the nonconservation of energy. The probability Wg % n ->s z ±i « thus becomes 
dominant for r=10 -5 and 10 -10 s, indicating that the magnitudes of the coefficients of Vsf 
are larger than those of the other coefficients. 

7. Conclusions 



Bearing in mind a single magnetic ion on a substrate or the ion trapped between a tip and 
a substrate, we derived a spin-atomic vibration interaction and a spin-flip Hamiltonian of a 
single atomic spin in a crystal field due to the surrounding ions. Here, the spin-atomic vibration 
interaction meant the interaction between the spin of the magnetic ion and the vibration of 
this ion. The spin-flip Hamiltonian represented terms containing spin-flip operators and no 
atomic vibration operators. 

In the derivation, we applied the perturbation theory to the following model: the crystal 
field potential energy in the equilibrium state was contained in the unperturbed Hamiltonian, 
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Fig. 14. Z component of the spin S z dependences of W^} n ^ Sz>n _ 1 /n, W§^ 2 n ^ Sz±ln _ 1 /n, 
Wl^ Sz ±i, n+ i/(n + 1), W§^ Sz>n _Jn, W§^ Sz , n+1 /(n + 1), W§l n _ Sz±hn of the 5=2 
system with huj=3\D\. Upper panel: r=10 -5 s. Lower panel: r=10~ 10 s. The energy levels of this 
system are shown in Fig. 13. 



while the spin-orbit interaction and the modulation of crystal field potential energy due to 
the vibration displacement were the perturbed terms. The model also took into account the 
difference in vibration displacement between the effective nucleus and the electrons, |ryAr n |, 
where Ar n is the displacement of the nucleus, and r\ (0 < r/ < 1) is a dimensionless quantity 
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Fig. 15. Z component of the spin Sz dependences of Wg^n-tSz n _i/n, W§ z 2 n ^ Sz±1 „_ 1 /n, 
Wl^ Sz ±i,n + i/(n + 1), Wl^s^/n, W|^ s , >n+1 /(n + 1), WjJ^s.ii,* of the 5=2 
system with hu=2\D\. Upper panel: t—10~ 5 s. Lower panel: t=1CP 10 s. The energy levels of this 
system are shown in Fig. 13. 



characterizing the difference. The unperturbed wave function ip^ (C = l " 5) was roughly set to 
be ip£= d,Q + Cd ^2m=i(^() dm + c p Yln=i Pn- Here, was the dominant 3d orbital, while d m 
and were the other 3d orbital and 4p orbital in the same magnetic ion, respectively. The 
coefficient for the d orbital, Cd, (the coefficient for p orbital, Cp,) represented the degree of d-d 
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hybridization (d-p hybridization). 

As results, the spin-atomic vibration interaction was obtained as Yli c i,iSi(~ a + 
a') + j{c2j,jSiSja + c^jjSjSja^), and the spin-flip Hamiltonian was expressed as 
Y^ij j c 4,i,jSiSj, with /, J=X, Y, Z. Here, Si was the / component of the spin 
S=(Sx, Sy, Sz), and a (a^) was the annihilation operator (creation operator) of the vibration 
of the magnetic ion. The coefficients of the spin-atomic vibration interaction c±j, C2j t j, and 
c 3,i,J contained 7], c p , Cd, and /, while that of the spin-flip Hamiltonian c^j t j had Cd, where / 
was a vibration frequency. 

The respective details about the spin-atomic vibration interaction, spin-flip Hamiltonian, 
and transition probability per unit time are written as follows: 

(1) Spin- atomic vibration interaction 

a. The spin-atomic vibration interaction appeared owing to the difference in vibration 
displacement and the d-p hybridization. Concretely speaking, this interaction was pro- 
portional to rj and also present for c p /0. The condition of c p ^0 was attributed to the 
result that the coefficients of this interaction contained the similar formula to the matrix 
element of the electric dipole transition. 

b. Most of the coefficients of the SiSja^ and SiSja^ terms (i.e., C2,/,j and csj y j) showed 
the 1/y/y dependence, indicating that these coefficients were proportional to the vi- 
bration displacement of the nucleus, Ar n (oc 1/y/J). In contrast, the coefficient of the 
Si(— a^ + ap term (i.e., c\j) had the v7 dependence because it was proportional to the 
momentum of the nucleus, Ap„ (oc y/J). Note that the coefficient of the conventional 
spin-phonon interaction was proportional to y/J, where this spin-phonon interaction 
meant an interaction between a single spin and the phonon of a lattice system, 

c. In the case of the Fe ion in the crystal field of tetragonal symmetry, the magnitudes 
of the coefficients of the spin-atomic vibration interaction became larger than that of 
the spin-phonon interaction by Mattuck and Strandberg in the specific region of /. The 
main reason for this was that the mass of the oscillator of the spin-atomic vibration 
interaction was smaller than that of the spin-phonon interaction. Here, the oscillator 
of the spin-atomic vibration interaction was the single magnetic ion, while that of the 
spin-phonon interaction consisted of all the ions in the unit cell. Note that the decrease 
in the mass of the oscillator increased the magnitude of vibration displacement, and 
furthermore the increase in the magnitude of the displacement enhanced the magnitude 
of the dominant term of the spin-atomic vibration interaction. 

(2) Spin-flip Hamiltonian 

a. The spin-flip Hamiltonian stemmed from only the d-d hybridization, and was indepen- 
dent of the atomic vibration. In other words, the spin-flip Hamiltonian existed only for 
Cdj^O, and did not contain / and rj. 
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b. In the case of the Fe ion in the crystal field of tetragonal symmetry, the magnitudes 
of the coefficients of the spin-flip Hamiltonian tended to be larger than those of the 
spin-atomic vibration interaction. This tendency was attributed to the result that the 
coefficients of the spin-flip Hamiltonian did not contain 77, whereas those of the spin- 
atomic vibration interaction were proportional to 77, where 77=0.05 was set in this paper. 

(3) Transition probability per unit time 

The transition probability per unit time was investigated for a simple model of the Fe 
ion. The unperturbed Hamiltonian was — \D\S Z + hwa^a, with 5=2, where D was an 
anisotropy constant, and oj (=27r/) was the angular frequency. As to —\D\S 2 Z , the en- 
ergy difference between the ground states and the first excited states was 3|D|, while 
that between the first excited states and the second excited state was \D\. The perturbed 
Hamiltonians were composed of the representative terms in the spin-atomic vibration in- 
teraction and spin-flip Hamiltonian. Note here that the magnitudes of the coefficients of 
the spin-atomic vibration interaction were set to be smaller than those of the spin-flip 
Hamiltonian (see 2b). The transition probability due to the spin-atomic vibration inter- 
action (e.g., SxSz(a + a T )) corresponded to the probability for the transition with the 
conservation of energy, depending on hio and the initial spin state (i.e., Sz)- In contrast, 
the transition probability due to the spin-flip Hamiltonian (e.g., SxSz) became the prob- 
ability for the transition with the nonconservation of energy. In the cases of hu)=3\D\ and 
specific Sz 's, the transition probability due to the spin-atomic vibration interaction took 
the largest value in all the probabilities because of the transition with the conservation 
of energy. In the case of Hoj=2\D\, all the transition probabilities were the probabilities 
for the transitions with the nonconservation of energy. The transition probability due to 
the spin-flip Hamiltonian thus became dominant, indicating that the magnitudes of the 
coefficients of the spin-flip Hamiltonian were larger than those of the other coefficients. 
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Appendix A: Difference in Displacement between Nucleus and Electron 

Let us consider a model in which a single electron is present in the a-spin shell (e.g., 
Fe 2+ : an electron exists in the down-spin shell). For this model, we evaluate the difference in 
vibration displacement between an effective nucleus and the electron in the cr-spin shell. The 
effective nucleus consists of a nucleus and core electrons, where the core electrons correspond 
to the electrons other than those in the cr-spin shell. We here use the theory of electronic 
polarizability in ref. 53. 
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Fig. A-l. Schematic illustration of the vibrating magnetic ion. The effective nucleus has the displace- 
ment Ar n in the x direction. The electron then moves in a circular orbit around the nucleus fixed 
at Ar n , and simultaneously experiences the crystalline electric field due to the surrounding ions. 
The x component of the position of the electron is written as Ar n + x. 

We consider a situation in which only the magnetic ion vibrates in the x direction, while 
the surrounding ions are rigid. Such a situation is schematically illustrated in Fig. A-l. In 
particular, this ion acts as a harmonic oscillator in the vicinity of a stable equilibrium position. 
This position corresponds to a local minimum point of a potential energy of the system. 
Here, the energy consists of "a Coulomb potential energy between the magnetic ion and the 
surrounding ions" and "a repulsive potential energy due to the Pauli exclusion principle and 
the overlap of electrons between the magnetic ion and the surrounding ions". 54 Now, the 
magnetic ion is regarded as a mass point with the mass of this ion, M. The position of the 
mass point corresponds to the gravity point of the magnetic ion, which may be, in other 
words, the position of the effective nucleus. When this mass point has the displacement from 
the equilibrium position, Ar n , the equation of motion of the mass point is written as 



where w n is the vibration frequency of the magnetic ion. 

When the effective nucleus is located at Ar n , the electron is considered to move in a 
circular orbit around "the nucleus fixed at Ar n " . In particular, we focus on the x component 
of the circular motion of the electron in the xy plane (see Fig. A-l). The x component of the 
position of the electron is written by Ar n + x. Here, x (— r c < x < r e ) is the x component 
of a vector that points from the nucleus to the electron, with r c being the electron radius. 
In addition, the electron experiences the crystalline electric field due to the surrounding ions. 
For this electron, we have the following equation of motion in the x direction: 



where u e is the angular velocity of the circular motion of the electron, — e (< 0) is the electronic 
charge, and Eion is the x component of the crystalline electric field. 

On the left-hand side of eq. (A-2), we note that |Ar n | (=|cj„Ar n |) is estimated to be much 



MAr n = -Mu^Ar, 



(A-l) 



m(Ar n + x) = —muj^x — eE\. 



(A-2) 
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smaller than \x\ (~ because of 1CT 12 < (uj n /uj c ) 2 < and 9.47 x 1CT 3 < |Ar n /x| < 

9.47 x 10 -1 . We here use w n =27r/ with 0.01 THz < / < 100 THz, and lo c of eq. (A-6). The 



magnitude of the displacement |Ar n | has been considered to be 9.47xl0~ 4 nm < y Ar 2 < 
9.47xl0~ 2 nm in Appendix . In addition, |x| has been roughly replaced with the square root 
of the time average of x 2 , i.e., r e /V2 (~0.1 nm), where r e corresponds to r n of eq. (A-9). 
Thus, by neglecting Ar n , eq. (A-2) is rewritten as 

mi = —muj 2 x — eE\ on . (A-3) 

Regarding Ei on , we consider a case in which the surrounding ions with the minus charge 
— Z- lon e (< 0) are respectively located at — R and R in the x axis as the point charges. Under 
the conditions of r e /i?<Cl, \x\/R<^.l, and \Ar n \/R <Cl, E- lon is obtained as 

-Z ion e j -(R-Ar n -x) R + Ar n + x \ 

Aon — — ; S — : : tttk + 



4vre [ [(R - Ar n - x) 2 + r 2 - x 2 } 3/2 [(R + Ar n + x) 2 + r% - .<• 



2 _ ™21 3 / 2 



»^(. + ^, (A4) 

where terms higher than the first order of r c /R, x/R, and Ar n /R have been neglected. Note 
that, in eq. (A-4), the effect of the repulsive potential energy due to the overlap of electrons 
has not been taken into account for simplicity. Substituting eq. (A-4) into eq. (A-3), we obtain 

( 2 , "^ion e2 \ -^ion e2 a I A 

rax = — muj e H — r x —^Ar n . (A-5) 



ireoR 3 J ireoR 3 

Here, by using the speed of electron, v=r e uj e , the expression for the equilibrium of force, 
mv 2 /r e =Z c ge 2 / (Ane^r 2 ) , and the Bohr's quantization condition r c mv=nh with n=l, 2, 3, • • ■, 
uj e is expressed as 

with a B =€oh 2 /(irme 2 ), where a B is the Bohr radius. 

The center position of the oscillation of the electron, x c , corresponds to the difference in 
vibration displacement between the nucleus and the electron. This x c is obtained by putting 
x=0 and x=x c in eq. (A-5): 



x c = -r)Ar n , (A-7) 

4(Z ioQ /Z cS ) (r n /Rf 
l + 4{Z io jZ c z){r n /Rf' 

r n = n 2 a B /Z cS , (A-9) 



with < rj < 1. Here, x c =0 is confirmed for Z[ on =0. 

In Fig. A-2, we show the Z[ on /Z e Q dependence of rj for R=2.5r n , 3r„, 4r n , and 5r n . The 
quantity r/ increases with increasing Z- lon /Z e R. In the case of Z- lon /Z e Q =0.2, we obtain 77=0.05 
for R=2.5r n , 77=0.03 for R=3r n , r/=0.01 for R=4r n , and 77=0.006 for R=5r n . 
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As an application, we consider r\ of the single Fe ion on the CuN surface. 10 The effective 
nuclear charge of the 3d orbital of Fe, Z e s, is obtained as Z e g=6.25 according to Slater's 
rules, 37 while Z- lon arising from the nearest-neighbor N is set to be Z- lon =lA on the basis of 
ref. 10. From Zi on /Z e ff=0.224, the largest rj is evaluated to be about 0.05. 




Appendix B: Magnitude of Displacement of Nucleus |Ar n | 

On the basis of eq. (A-l), we evaluate the magnitude of vibration displacement of the 
effective nucleus, |Ar n |. Here, |Ar n | is defined as the square root of the time average of Ar 2 . 
The solution of eq. (A-l) is first obtained as 

Ar n = Acos(w n t) + Bsm(uj n t), (B-l) 

where A and B are coefficients. The total energy is then written as 

E = iMAf* + \mu\AvI 

= l -Mu>l{A 2 + B 2 ) 

= Mc^ATf. (B-2) 
Here, Ar^ is the time average of Ar 2 , given by 

ATI = lim \ f Arl&t = \{A 2 + B 2 ). (B-3) 
t'-s>oo v Jq 2 

On the other hand, the energy of a quantum harmonic oscillator is expressed as ((n) + 
l/2)huj n . Here, (n) is the thermal average of the vibrational quantum number, given by (n) 
= [exp(^ n //c B r) — l] -1 , where T is the temperature and fee is the Boltzmann constant. 
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We now equate E of eq. (B-2) with the above energy, that is, 

1 



Mu 2 n Arl =((«) + -) Hun. (B-4) 



From eq. (B-4), we can obtain \J Ar r 2 r At T=0 K, \J Ar 2 is written as 



The quantity y^ Ar 2 decreases with increasing M and w n . 

We also rewrite eq. (A-l) as MAf n = — CAr n , with C = Mw„. Here, C is a force constant 
that corresponds to the second derivative of the interatomic potential with respect to the 
displacement of the magnetic ion. 56 Substituting uj q =^C/M into eq. (B-5), we obtain 



The quantity y Ar 2 decreases with increasing M and C. 

As an application, we consider \J Ar 2 at T=0 K for Fe 2+ , where M of Fe 2+ is 
M«(26+30)x(1.67xl(T 27 ) kg. By using eq. (B-5), yffir* is evaluated to be 9.47xl0~ 4 nm 
< y^Arf <9.47xl0~ 2 nm for 0.01 THz < / < 100 THz, where co n = 2tt/. 

Appendix C: Expression of Coefficients for Fe Ion 

We describe the expressions of coefficients of the spin-atomic vibration interaction and 
the spin-flip Hamiltonian of the Fe ion (Fe 2+ ) in the crystal field of tetragonal symmetry. For 
convenience, the coefficients are represented by the following dimensionless quantities: A^/A, 
A^/A, G C> F/D, A'j j.^/D, K% v JD, A^JD, and Afj/D. 

We here focus on one down-spin electron as described in §4. The energy levels are assumed 
to be Ej of Fig. 2. In addition, we use the notation of \^i)=\xy), \^2)=\yz), \^^)=\xz), 
|* 4 )=|x 2 -y 2 ), |* 5 )=|3z 2 -r 2 ), |di)=|xy), \d 2 )=\xz), \d 3 )=\yz), \d 4 )=\x 2 -y 2 ) , \d 5 )=\3z 2 -r 2 ) , 
|pi)=|x), \p2)=\y), and \ps)=\z). As to Ej, ej, cf^, Cp}, and Cj, we have the notation and 
parameter setup in §4.2. 

In the below-mentioned calculations, we will use the following matrix elements: 

(xz\L x \xy) = i, (C-l) 

(x 2 - y 2 \L x \yz) = i, (C-2) 

{3z 2 - r 2 \L x \yz) = V3i, (C-3) 

(yz\L y \xy) = -i, (C-4) 

(x 2 — y 2 \L y \xz) = i, (C-5) 

(3z 2 - r 2 \L y \xz) = ->/3i, (C-6) 



43/55 



J. Phys. Soc. Jpn. 





Pttt t Pa ppt? 


/ r 2 _ ,/ 2 | r | T7 A — _9i 




ll'? /, f/TM — — 1 

\jy6\ij z \y6] — 1, 




(z\L x \y) = i, 


(C-9) 


(x\L y \z) = i, 


(C-10) 


(y\L z \x) = i, 


(C-ll) 



with (dj/\L\dj)= —(dj\L\dj/) and (j)j/|L|j)j)= — (pj\L\pj>) 
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C. 1 Principal axis coordinate system 

When an electron in the crystal field of tetragonal symmetry has Ej in Fig. 2, we find 
that X, Y, and Z are written as X=x, Y=y, and Z=z, respectively. 

(0) 



Using eqs. (C-l) - (C-8), we obtain C^l of eq. (90) as 





/ 


AO) 

'~'X,X 





\ 


£(0) = 







£ (0) 
l ~"y,y u 











£g] 


/ 



where 



4°j» = ^ 



(xy|L x |x2)(x2|.L x |a;y) 



^4 

A,' 



V>V o - t, A ' 



&xy &yz 



(xy\L z \x 2 - y 2 ){x 2 - y 2 \L z \xy) 



C= (l + 4cg + 3<£ 



G-xy e x 2 —y 2 
2\-V2 



\C A 



with C^ u =0 for \i^v. Here, we have roughly set e\ — ej=E\ — Ej, taking into account | Aei — 
Aej|<C|ei — ej\ for j ^ 1 (see eq. (45)). 

Since is a diagonal matrix, we obtain 0x,x=&Y,y=9z,z = Q m e Q s - (94) and (95), and 
the other 9i^s are ir/2. Namely, X, Y, and Z become X=x, Y=y, and Z=z, respectively. 



(C-12) 



(C-13) 

(C-14) 

(C-15) 
(C-16) 



C.2 D, E, and F 

Using eqs. (103) - (105) and (C-13) - (C-15), we obtain D, E, and F as 

o = x 2 ( c {0) - flk±f£l 1 - ^91 ( 



z,z 



AO) _ AO) 
E = X 2 Lx ' X 2 Ly ' Y =0, 

AP) , £ (o) 



A"' 



(C-17) 

(C-18) 
(C-19) 
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with x =C x °}x, /2yy=/2^y, and 4°z = ^°i' The bistability (i.e., D < 0) is obtained for 
A m /A s < 4. The coefficient D comes close to — 3A 2 C 4 /A m as A m /A s comes to 1. In addition, 
we find that the conditions of A s ~ A m and a small A m as well as a large |A| enhance \D\. 
The coefficient E becomes owing to x =Cy\- 



C.3 A«/A 

Using eqs. (55), (7), and |^i)=|xy), we have 



A 



with fi, £=x, y, z. By using \xy) in eq. (122), (xy\/i\xy) is obtained as follows: 



(C-20) 



(xy\x\xy) = C 2 c p (^(xy\x\y) + (y\x\xy)^j + c d ^(xzjxjz) + (z|x|xz) + (x 2 — y 2 \x\x) + (x|x|x 2 — y 2 ) 



+ {3z 2 — r 2 \x\x) + (x\x\3z 2 — r 2 )^j 



(C-21) 



(xy\y\xy) = C 2 c p ({xy\y\x} + (x\y\xy)j + c d ((yz|y|z) + (z\y\yz) + (x 2 - y 2 |y|y) + {y\y\x 2 - y 2 ) 



+ (3z 2 - r 2 \y\y) + (y\y\3z 2 - r 2 )) 



(C-22) 



(xy|z|xy) = C 2 CdCp {{xz\z\x) + (x|.z|xz) + (yz|z|y) + (y\z\yz) + (3z 2 - r 2 \z\z) + (z|z|3z 2 - r 2 )^ , 

(C-23) 

where (xy\/i\xy) vanishes at c p =0. The matrix element (pji\(j,\dj) is obtained using eqs. (125) 
- (132): 

(x\y\xy) = -8.31 x 10~ 12 , (C-24) 
(x|x|3z 2 -r 2 ) =4.80 x 10~ 12 , (C-25) 
(z\z\3z 2 - r 2 ) = -9.59 x 1(T 12 , (C-26) 
(x|x|x 2 -y 2 ) = -8.31 x 10~ 12 , (C-27) 
(y\y\x 2 -y 2 ) = 8.31 x 10~ 12 , (C-28) 

where (dj\fi\pj/)=(pji\fi\dj) . 



c.4 Ayx 

Using eqs. (80) and (7), we obtain Aj^/A as 



A 



(C-29) 



'M/iV 2 

for |Lt, £=x, y, z. The matrix element (xy\L\^j) is expressed by 

(xy\L x \yz) = C 2 c d (xy\L x \xz) + c d ((x 2 - y 2 \L x \yz) + (3z 2 - r 2 \L x \yzfj + c 2 {{y\L x \z) + (z|L x |y)) 
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+c 2 d {(xz\L x \xy) + (yz\L x \x 2 - y 2 ) + (yz\L x \3z 2 - r 2 ))] = s/3c d (l - c d )C 2 i, (C-30) 



(xy\L x \xz) = C 



{xy\L x \xz) + c 2 {{y\L x \z) + {z\L x \yfj + c 2 d {{xz\L x \xy) + {x 2 - y 2 \L x \yz) 
+W\L x \x 2 - y 2 ) + (3z 2 - r 2 \L x \y~z~) + W\L x \3z 2 - r 2 ))] = (c 2 d - 1)C% (C-31) 
(xy|L x |x 2 - y 2 ) = C 2 c d ((xy\L x \xz) + (yz\L x \x 2 - y 2 )) + c 2 p {{y\L x \z) + (z|L x |y)) 

+c 2 (^z~\L x \x~y) + (x 2 -^^) + (S^-r 2 !^^) + W\L x \3z 2 - r 2 )) 

= 2 C(J (c d - l)C 2 i, (C-32) 
(xylL^z 2 - r 2 ) = C 2 [c d ({xy\L x \x~£j + M^|3^ 2 - r 2 )) + c 2 ( (y] L x ^) + (4^h/>) 

+ c d ((^l-L^x?/) + {x 2 - y 2 \L x \yz) + (y^L^x 2 - y 2 ) + (3z 2 - r 2 \L x \yz)^j 

= (1 + ^)c d ( C(i - l)C 2 i, (C-33) 
(xy|L y |yz) = C 2 (xy\L y \yz) + c 2 ((x|L y |z) + (z\L y \x)^j + c 2 , ((yzlL^xy) + (x 2 - y 2 \L y \xz) 

+Jx~z~\L y \x 2 - y 2 ) + (3z 2 -r 2 \L y \x~z) + (^L y |3z 2 - r 2 )) = (1 - c^)C 2 i, (C-34) 
c d ((xy\L y \yz) + (x 2 - y 2 \L y \xz) + (3z 2 - r 2 |L y |x2;)^ + c 2 ((xl-L^z) + 



(xy\L y \xz) = C 



+c 2 d {[yz\L y \xy) + {xz\L y \x 2 - y 2 ) + (xz\L y \3z 2 - r 2 ))] = (2 - >/3)cj(l - c d )C\ (C-35) 



(xylL^x 2 - y 2 ) = C 2 c d ((xy|L y |yz) + (xzjL^lx 2 - y 2 )) + c 2 ((x|L y |;<:) + (z|L y |x)) 

+C 2 , [{yz\L y \xy) + (x 2 - y 2 \L y \xz) + (3z 2 - r 2 |L y |xz) + (xz\L y \3z 2 - r 2 )^j 



= 0, 
(C-36) 



(xy\L y \3z 2 - r 



C 2 



c d \ {xy\Ly\yz) + {xz\L y \3z 2 - r 2 }) + c 2 ^(xl-L^z) + (^L^x) 

+ c d {{vA L y\ x y) + - V 2 \ L y\ xz ) + (xz|L y |x 2 - y 2 ) + (3z 2 - r 2 \L y \xzf) 
= (1 + v/3> d (l - c d )C 2 i, 

((xy|L 2 |x 2 -y 2 ) + (^i|L z |yz)) + c 2 (jx\L g \y) + (feR) 

+ c d ((^ 2 - y 2 l^ky) + (y^l^l^)) = c d (i - c d )c\ 

((xy\L z \x 2 - y 2 ) + ~{y~z]L z \xz)^ + c 2 ((4^b> + W\ L z]x}) 



(xy\L z \yz) = C 



(xy\L z \xz) = C 

+c 2 (jxV\L z \y^ + (x 2 - y 2 \L z \x~y))] = 3c d (l - c d )C 2 \, 
(xy\L z \x 2 - y 2 ) = C 2 [(xy|L 2 |x 2 - y 2 ) + c 2 (jx\L x \y) + Jy\L z \x) ) 

+c 2 d ((x 2 -y 2 |L 2 ^) + M^.P) + (Ri z R)] = 2(1 - c\)C\ 
{xy\L z \3z 2 - r 2 ) = C 2 [c d (xy\L z \x 2 - y 2 } + c 2 p (jx~\L z \y) + ~&\L z \xf) 



(C-37) 



(C-38) 



(C-39) 



(C-40) 
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+c 2 (<x 2 -y 2 |L^> + JxY\L z \y~^ + <^M)] = 2c d (l - c d )C 2 i, (C-41) 

where (*j|L|*j)=0 and (^ / |L|4' i )=-(4' i |L|4' i /) for j/j'. 42 The matrix element (*j|£|xy) 
is obtained as 

(xz|x|xy) = C 2 c p ^(xz|x|z) + {y\x\xy)^j + ^(z|x|xz) + (xy\x\y) + (x|x|x 2 — y 2 ) 

+ {x 2 — y 2 \x\x) + (x|x|3z 2 — r 2 ) + (3z 2 — r 2 |x|x)^ , 
(xz\y\xy) = C 2 c p (x\y\xy) + c d {{xy\y\x) + (z\y\yz) + (yz\y\z) + (y|y|x 2 - y 2 ) + (x 2 - y 2 \y\y) 



(C-42) 



+ (y|y|3z 2 - r 2 ) + (3z 2 - r 2 |?/|y)) 
(xz|z|xy) = C 2 c p (xz\z\x) + (^(x\z\xz) + + (yzlzjy) + (z\z\3z 2 — r 2 ) + (3z 2 — r 2 \z\z)^j 

(yz\x\xy) = C 2 c p (y\x\xy) + c d (^(xy\x\y) + (z\x\xz) + (xz|x|z) + (x|x|x 2 — y 2 ) + (x 2 — y 2 |x|x) 

+ (x|x|3z 2 — r 2 ) + (3z 2 — r 2 |x|x)^ , 
(yz\y\xy) = C 2 c p ((yz\y\z) + (x\y\xyf) + c d {{z\y\yz) + (xy\y\x) + {y\y\x 2 - y 2 ) + {x 2 - y 2 \y\y) 

+ {y\y\3z 2 - r 2 ) + (3z 2 - r 2 \y\y)^j , 
(yz\z\xy) = C 2 c p (yz\z\y) + c d (^(y\z\yz) + (xz\z\x) + (x\z\xz) + (z|z|3z 2 — r 2 ) + (3z 2 — r 2 \z\z)^j 



(C-43) 



(C-44) 



(C-45) 



(C-46) 



(C-47) 



(x — y |x|xy) = C c p i(x — y \x\x) + {y\x\xy) J + c d ( (x|x|x 2 — y 2 ) + (xy\x\y) + (xz|x|z) + (z|x|xz) 



+ (x|x|3z 2 — r 2 ) + (3z 2 — r 2 |x|x)^ 



(C-48) 



(x -y\y\xy) = Cc p ( (x - y \y\y) + (x\y\xy) ) + c d ( (y\y\x 2 - y 2 ) + (xy\y\x) + (z\y\yz) + (yz\y\z 



+ (y\y\3z 2 - r 2 ) + (3z 2 - r 2 |y|y)) 



(C-49) 



(x 2 — y 2 \z\xy) = C 2 c d c p (^{x\z\xz) + {xz\z\x) + (y\z\yz) + (yz\z\y) + (z|z|3z 2 — r 2 ) + (3z 2 — r 2 |z|z)^ , 



(3z 2 — r 2 |x|xy) = C 2 c p (^(3z 2 — r 2 |x|x) + (y\x\xy)^j + c d ^(x|x|3z 2 — r 2 ) + (xy\x\y) + (xz\x\z) 



+{z\x\xz) + (x 2 — y 2 |x|x) + (x|x|x 2 — y 2 ) 



(3z 2 - r 2 \y\xy) = C 2 c p {{3z 2 - r 2 \y\y) + (x\y\xy)^ + c d (^(y\y\3z 2 - r 2 ) + {xy\y\x) + {yz\y\z) 



+ {z\y\yz) + (x 2 - y 2 \y\y) + (y|y|x 2 - y T 



(C-50) 



(C-51) 



(C-52) 



(3z 2 — r 2 \z\xy) = C 2 c p {3z 2 — r 2 \z\z) + c d (^(z\z\3z 2 — r 2 ) + (xz\z\x) + (x\z\xz) + (yz\z\y) + 
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where (* J -,|e|* i )=(*il^l*i')- Here, (pf^dj) is given by eqs. (C-24) - (C-28). 



(C-53) 



C.5 Afj/D 

Using eqs. (106) and (103), we obtain Aj F j/D as 



ASF 



D ~ r (0) - cr (0) 4- r (0) w?' 

with /, J=X, Y, Z. Here, cf z - {C% ] x + 4? Y )/2 (see eq. (C-17)) becomes 



(C-54) 



r (o) , r (o) 

^z,z 



1 



C 4 . 



(C-55) 



2 \ A s A m , 

Using eqs. (C-30) - (C-41), we obtain C'j j of eqs. (97) (or C! of eq. (91)) as 

(xy\L x \yz)(yz\L z \xy) (xy\L x \xz)(xz\L z \xy) (xy\L x \x 2 - y 2 )(x 2 - y 2 \L z \xy) 



£'x,Z — £>'x,z 



E — E 

'-'xy J -'yz 



+ 



E — E 
1 - J xy '-'xz 



+ 



E xy E x 2_ y 2 



{xy\L x \3z 2 - r 2 ){3z 2 - r 2 \L z \xy) 



E xy — E 3z 2_ 



_£(0) 



C 4 c d (c d - 1) 



(3 - y^c 2 + V3c d - 3 | 4(c 2 -l) | 2(1 + ^3)^-1) 



(C-56) 



, , _ (xy\L y \yz)(yz\L z \xy) (xy\L y \xz){xz\L z \xy) (xy\L y \3z 2 - r 2 )(3z 2 - r 2 \L z \xy) (0) 

^Y,Z — L 'y,z — ^ ^ 1 ^ ^ 1 ~ ~ ■> 



E xy E yz 



E — E 

'-'xy '-'xz 



E X y - E 3z 2_ 



V--Z 



C 4 c d (c d - 1) 



(-7 + 3V3)c 2 + (6 - 3y/3)c d + 1 -(2 + 2y/ 3)cg + (2 + 2>/3)cd 



A, 



A, 

, _ , _ (xy\L x \yz)(yz\L y \xy) (xy\L x \xz)(xz\L y \xy) (xy\L x \3z 2 - r 2 )(3z 2 - r 2 \L y \xy) _ (0) 

*~X,Y ~ >->x,y ~ jp — tp + f _ jr + B7 _ B!_ „ „ ^.1/ 



(C-57) 



-X,Y ^X,y TP TP T JP — JP ^ JP. — En 9 9 

'-'xy J -'yz '-'xy '-'xz '-'xy '-'J,z i —r i 



-'xy '-'xz 



C 4 c d (c d - 1) 



(2-2^/3)^ + 2^/3-2 (l + V3) 2 (cg-c d ) 

A, + A; 



(C-58) 



E x y E xz 



f., ~, ___ (xy\L x \yz)(yz\L x \xy) (xy\L x \xz)(xz\L x \xy) (xy\L x \x 2 - y 2 )(x 2 - y 2 \L x \xy) 

x ' x x ' x E — E E — E E - E o o 

'-'xy ^yz '-'xy '-'xz '-'xy Lj x 1 —y 1 

, {xy\L x \3z 2 - r 2 ){3z 2 - r 2 \L x \xy) 



E X y — E 3z 2_ 



_£(0) 



C A c d 



Ac d + 2 4c d (1 + V3) 2 c d 



A, 



A r . 



A, 



+ (cj - 2) 



4c 2 + 2c d + 1 _ 4c| _ (1 + y/3) 2 c; 
A™ 



A, 



A, 



(C-59) 



, (xy|L y |y«)(y«|L y |xy) (xy|L ?/ |xz)(xz|L ?/ |xy) (:n/|Lj,|3z 2 - r 2 )(3z 2 - r 2 \L y \xy) (0) 

— L 'y,y — ^ ^ 1 ^ ^ 1 ^ ^ 7/ 



2_ r 2\^g^2 „2| 



jp. — E ' E — En 9 - ~" y,y 

'-'xy '-'xz '-'xy '-'Zz 1 - 

, (8-4^3)c d + 2 (l + v 7 ^) 2 ^ 



A, 
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(8 - AVS)c 2 d + 2c d + 1 (1 + v/3) 



2„2 



(xy\L z \yz)(yz\L z \xy) (xy\L z \xz)(xz\L z \xy) (xy\L z \x 2 - y 2 )(x 2 - y 2 \L z \xy) 



(C-60) 



B — B 



+ 



E — E 



+ 



E X y E x 2_y2 



{xy\L z \3z 2 - r 2 ){3z 2 - r 2 \L z \xy) _ (0) 

I? IT 1 Z,Z 

^xy ~ ^3z 2 -r 2 



C 4 c d 



10c rf 4(c d + 2) 4c d 

a7"^ A7 + (cd " 2) 



lOj 4fa + l) 2 4c 2 
A s A m A t 



(C-61) 



with >Cx°z=>Cj ( °2=jCi°y=0. In addition, we have C' ZX =C' ZX =C' XZ , C! ZY =C! zy =CJ yz 
CJ Y x =C! y x =C xy . Note that Cl 1 j becomes at c d =0. 

C.6 G e and F jD 

The coefficient is given by eq. (110): 



G 5 = 2t7 



<3^ 



ft 



(C-62) 



(r- 3 ) V 2Mw ? ' 

with £=x, y, z. Here, (r -3 ) and (3£r~ 5 ) are regarded as the respective ground-state expectation 
values: 

( r - 3 ) = ( xy \r^\xy) 

= C 2 {xy\r~ 3 \xy) + c 2 ((^r _3 ]x) + Jy\r~ 3 ]y) + (i|r~ 3 |i)^) + c 2 , (jyz\r- 3 ]yz) + (xi|r~ 3 |xi) 
+(x 2 -y 2 |r~ 3 |x 2 -y 2 > + (3z 2 -r 2 |r" 3 |3z 2 - r 2 )) ] , (C-63) 



(3xr 5 ) = (xy|3xr 5 |xy) 



(xy|3xr + (y|3xr °\xy)) + c d [(x 2 — y 2 \3xr °|x) + (x|3xr 5 |x 2 — y 2 ) 



+ (3z 2 — r 2 |3xr 5 |x) + (x|3xr 5 |3z 2 — r 2 ) + (xz|3xr 5 \z) + (z\3xr 5 \xz)^j 
(3yr~ 5 ) = (xy\3yr~ 5 \xy) 



(C-64) 



C Cr) 



{{xy\3yr 5 |x) + (x|3yr 5 |xy)) + c d ((x 2 - y 2 \3yr 5 \y) + {y\3yr 5 |x 2 - y 2 ) 



(C-65) 



+ (3z 2 - r 2 |3yr~ 5 |y) + (y|3yr~ 5 |32: 2 - r 2 ) + (y2;|3yr _5 |2;) + (z|3yr~ 5 |?/2:)^ 
(3zr~ 5 ) = (xy\3zr^\xy) 

= C 2 c d c p (^(xz\3zr~ 5 \x) + (x\3zr~ 5 \xz) + (yz\3zr~ 5 \y) + (y\3zr~ 5 \yz) 

+ (3z 2 - r 2 \3zr^\z^ + Jz~\3zr~ 5 \3z 2 - r 2 )) , (C-66) 



with 



(x|r" 3 |x) = 2.42 x 10 30 , 
{xy\r~ 3 \xy} = 4.08 x 10 30 , 



(C-67) 
(C-68) 
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{3z 2 - r 2 \r~ 3 \3z 2 - r 2 ) = 4.06 x 10 30 , (C-69) 

(x 2 - y 2 \r- 3 \x 2 - y 2 ) = 4.06 x 10 30 , (C-70) 

(xy\3xr- 5 \y) = 1.26 x 10 41 , (C-71) 

(x 2 - y 2 \3xr- 5 \x) = 1.26 x 10 41 , (C-72) 

(x 2 - y 2 \3yr- 5 \y) = -1.26 x 10 41 , (C-73) 

{3z 2 - r 2 \3xr- 5 \x) = -7.26 x 10 40 , (C-74) 

{3z 2 - r 2 \3zr~ 5 \z) = 1.45 x 10 41 , (C-75) 

where eqs. (C-67) - (C-75) have been obtained by using eqs. (125) - (132). 



C.7 A' IfJ . 4 /D 

Using eqs. (Ill) and (103), we have 



h 



i.j 



(C-76) 



D - (r-3) \l2M^ C f z _ {C f x + C W Y)/2 > 

for I, J=X, Y, Z and £=x, y, z. This Aj j ^/D is also represented by A' IJ JD=G^A^ J /D 
(see eqs. (110) and (C-54)). 



Using eqs. (C-30) - (C-41), we express A^. u ^/D of eq. (81) with fx, v, £=x, y, z as 



D 



A ± 
D 



D 



with 



VSc d (l - c d )\(yz\v\xy) ± (c 2 d - l)i(xz\v\xy) ± 2c d (c d - l)i(x 2 - y 2 \v\xy) 



-Ac 



(l + V3)c d (c d -l)i(3z 2 -r 2 |i/M ± 



± 



A, A, 



(C-77) 



(1 - <j)i(j/z|i/|ay) ± | (2 - \/3)cd(l - Cd)i(az|i/|a;y) ± 

=ai k ^ + =ai K ^ 



(1 + V3)c d (l - c d )i(3z 2 -r 2 \u\xy) ± 



1 4 )C 2 
aT" a7 ,c ' 



(C-78) 



c d (l -Cd)i(^|i/|a;y) ± 3c d (l - c d )i(az|i/|xy) ± 2(1 — cg)i(x 2 — j/ 2 |z^|xj/) ± 
-A s -A s -A m 



, 2c d (l - c d )i(3z 2 -r 2 |^|x|/) ± 



1 



c 2 , 



(C-79) 



_i_ i ( m Mfvjjt 



' 7i m 
2Mu{h 



A., 



(C-80) 
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for u=s, m, and I. The matrix element 1^-) is given by eqs. (C-42) - (C-53). We also 

note that satisfies Aj^-A^ = (Aj^)* for M, J/, z 44 where Aj^ is eq. 

(84). 
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Table I. We show the presence or absence of the spin-atomic vibration interactions ^g Q]j) 

u soO,sol,i' w S),bo1,2' and ^00,802' and the spin-hip Hamiltonian v^ sqQ 2 for each set of V , c p , 
Cd- Here, ei — ej=E\ — Ej is set by taking into account |Aei — Ae.,|<sj;|ei — ej\ for j ^ 1. The 
spin-atomic vibration interaction (spin- flip Hamiltonian) is named SA (SF). The presence and 
absence are represented by y/ and x , respectively. The operators and coefficients are also written 
for reference. The coefficient A^.^ contains both terms with l/y/fe and those with \fJl- The 
coefficient Af^ does not contain 77, /j, and M. 



Name SA SA SA SA SA SF 

Enerev V {1) V {2) v {2) v {2) V {2) v {2) 

energy v &o2 soO,k v soO,sol,l v so0,sol,2 *so0,so2 v so0,so0,2 

eq. (116) eq. (114) eq. (108) eq. (109) eq. (112) eq. (102) 

Operator Si(—a^ + a^) Si(—a^ + a^) S](a^ + a^) SiSj(a^ + a^) SjSja^ SjSj 

SiSja\ 

Coefficient A?* A* GtD A' A± A? F , 

G^E 
G(.F 

x [h x [k V>? V>? 7]X 2 2 

V m V m yfWl \[Wi Vm 

eq. (55) eq. (80) eq. (110) eq. (Ill) eq. (81) eq. (106) 



77=0 

c d =0, c p =0 x 

c rf =0, Cp/0 x 

Cd^O, c p =0 x 

Cd/0, Cp/0 x 



x 

X 
X 
X 



X 
X 
X 
X 



X 
X 
X 
X 



X 
X 
X 
X 



V 
V 

















c d =0, c p =0 


X 


X 


X 


X 


X 


X 


Q=0, Cp/0 




V 


V 


X 




X 


c p =0 


X 


X 


X 


X 


X 




Cp/0 








V 


V 





52/55 



J. Phys. Soc. Jpn. Full Paper 

References 

J. Villain, F. Hartmann-Boutron, R. Sessoli, and A. Rcttori: Europhys. Lett. 27 (1994) 159. 
M. Misiorny and J. Barnas: Phys. Rev. B 77 (2008) 172414. 

A. Ardavan, O. Rival, J. J. L. Morton, S. J. Blundell, A. M. Tyryshkin, G. A. Timco, and R. E. 
P. Winpenny: Phys. Rev. Lett. 98 (2007) 057201. 

A. Caneschi, D. Gatteschi, R. Sessoli, A.-L. Barra, L.-C. Brunei, and M. Guillot: J. Am. Chem. 
Soc. 113 (1991) 5873. 

R. Sessoli, H.-L. Tsai, A. R. Schake, S. Wang, J. B. Vincent, K. Folting, D. Gatteschi, G. Christou, 
and D. N. Hcndrickson: J. Am. Chem. Soc. 115 (1993) 1804. 

P. Politi, A. Rettori, F. Hartmann-Boutron, and J. Villain: Phys. Rev. Lett. 75 (1995) 537. 
L. Thomas, F. Lionti, R. Ballou, D. Gatteschi, R. Sessoli, and B. Barbara: Nature 383 (1996) 145. 
K. Takcda and K. Awaga: Phys. Rev. B 56 (1997) 14560; K. Takcda, K. Awaga, and T. Inabe, 
Phys. Rev. B 57 (1998) R11062. 

A.-L. Barra, P. Debrunner, D. Gatteschi, Ch. E. Schulz, and R. Sessoli: Europhys. Lett. 35 (1996) 
133. 

C. F. Hirjibehedin, C.-Y. Lin, A. F. Otte, M. Ternes, C. P. Lutz, B. A. Jones, and A. J. Heinrich: 
Science 317 (2007) 1199. 

C. F. Hirjibehedin, C. P. Lutz, and A. J. Heinrich: Science 312 (2006) 1021. 
J. Park, A. N. Pasupathy, J. I. Goldsmith, C. Chang, Y. Yaish, J. R. Pctta, M. Rinkoski, J. P. 
Scthna, H. D. Abruha, P. L. McEuen, and D. C. Ralph: Nature 417 (2002) 722. 

D. Rugar, R. Budakian, H. J. Mamin, and B. W. Chui: Nature 430 (2004) 329. 
N. Tsukahara, K. Noto, M. Ohara, S. Shiraki, N. Takagi, Y. Takata, J. Miyawaki, M. Taguchi, A. 
Chainani, S. Shin, and M. Kawai: Phys. Rev. Lett. 102 (2009) 167203. 

Y. Yayon, V. W. Brar, L. Senapati, S. C. Erwin, and M. F. Crommie: Phys. Rev. Lett. 99 (2007) 
067202. 

C. Etz, J. Zabloudil, P. Weinberger, and E. Y. Vcdmcdenko: Phys. Rev. B 77 (2008) 184425. 

D. Serrate, P. Ferriani, Y. Yoshida, S.-W. Hla, M. Menzel, K. von Bcrgmann, S. Heinze, A. Kubetzka, 
and R. Wiesendanger: Nat. Nanotechnol. 5 (2010) 350. 

S. Kokado, K. Ueda, K. Harigaya, and A. Sakuma: Phys. Rev. B 76 (2007) 054451. 
S. Kokado, K. Harigaya, and A. Sakuma: Phys. Status Solidi C 6 (2009) 2113. 
For example, see S. Amasha, K. MacLean, I. P. Radu, D. M. Zumbu, M. A. Kastner, M. P. Hanson, 
and A. C. Gossard: Phys. Rev. Lett. 100 (2008) 046803. 

A. Zhao, Q. Li, L. Chen, H. Xiang, W. Wang, S. Pan, B. Wang, X. Xiao, J. Yang, J. G. Hou, and 
Q. Zhu: Science 309 (2005) 1542. 

R. Z. Huang, V. S. Stepanyuk, A. L. Klavsyuk, W. Hergert, P. Bruno, and J. Kirschner: Phys. Rev. 
B 73 (2006) 153404. 

J. H. Van Vlcck: Phys. Rev. 57 (1940) 426. 
J. H. Van Vleck: J. Chem. Phys. 7 (1939) 72. 

R. D. Mattuck and M. W. P. Strandberg: Phys. Rev. 119 (1960) 1204. 
R. Orbach: Proc. Phys. Soc. 77 (1961) 821. 
A. M. Stoneham: Proc. Phys. Soc. 85 (1965) 107. 

E. B. Tucker: Proc. IEEE 53 (1965) 1547. 



53/55 



J. Phys. Soc. Jpn. 



Full Paper 



29) A. M. Stoneham: Theory of Defects in Solids (Oxford University Press, Oxford, 1975) Chap. 14. 
In particular, see Table 14.1. 

30) F. Hartmann-Boutron, P. Politi, and J. Villain: Int. J. Mod. Phys. B 10 (1996) 2577. 

31) D. Gatteschi, R. Sessoli, and J. Villain: Molecular Nanomagnets (Oxford University Press, Oxford, 
2006) Chap. 5. 

32) In the spin-phonon interaction, 25 ' 30 each magnetic ion (or each magnetic molecule) in the crystal 
has a single spin, where the exchange interactions between the spins are neglected. The spin- 
phonon interaction originally represents an interaction between the spin located at the origin and 
the phonon of a lattice system. 

33) K. Liu and S. Gao: Phys. Rev. B 74 (2006) 195433. 

34) On the basis of the Coulomb energy between the 3d electrons, Vcouiomb = J2e,>£' e ' V (^ 7Te o\ r £~ r £' \), 
we show that the electrons prefer parallel spin alignment in the case of the weak crystal field. 
We now consider a case in which two electrons occupy the (C,cr) and (C',c') states, respec- 
tively, where ( (cr) represents the ( orbital (a spin). With the use of the multiorbital Hub- 
bard model, the energy may be roughly written as 2ed ( + U^^ for ( = Q' and a ^ a', and 
£d ( + £d c , + £^c,C — Jc,C f° r C C- Here, £<j { (e<i { ,) is an energy of the £ orbital ((' orbital), 
where sa ( < , (see eq. (39)). In addition, U^g is an intra-atomic Coulomb energy defined by 
U^x'= J J d^(£)d^,(£')e 2 / (4:ire \rt — re'\)d^(£)d^(£')dredrt' (> 0), and J^/ is an intra-atomic ex- 
change energy defined by Jc£>=f J d* i -{£)d*^{£')e 2 / {Ai:e \ri - re'\)d^(£')d^' (£)drtdre> (> 0), where 
d^(£) represents the 3d orbital (see eq. (36)). We assume the cases of [7^ = U^^i and the weak 
crystal field with Sd c , — £d ( < J(,('' where — Ed c is the crystal field splitting energy. We then 
obtain the relation of Ed ( + £d ( , + f^c.C — ^C,C < ^ £d c + ^C,C> which reduces to e<j , — < J(X'- 
The electrons thus occupy the orbitals, keeping their spins parallel. For the multiorbital Hubbard 
model, see L. M. Roth: Phys. Rev. 149 (1966) 306. 

35) K. Yosida: Theory of Magnetism (Springer Series, New York, 1998) Chap. 1. 

36) According to Slater's rule, 37 Z c s of Fe is evaluated to be 6.25. 

37) J. C. Slater: Phys. Rev. 36 (1930) 57. 

38) For example, see Appendix . As found from eq. (B-6), the magnitude of vibration displacement 
decreases with increasing mass and force constant. When the vibration of the surrounding ions is 
considered, the magnetic ion (the surrounding ions) in this Appendix should be replaced with the 
surrounding ions (the other ions adjacent to them in the substrate). 

39) M.-C. Desjonqueres and D. Spanjaard: Concepts in Surface Physics (Springer- Verlag, Berlin, 1995) 
2nd cd., p. 505. 

40) See chap. 3 in rcf. 35. 

41) M. H. L. Prycc: Proc. Phys. Soc. A 63 (1950) 25. 

42) We here use the wave functions of 3d orbitals, i.e., c£i =0^/3,2 (?*), d2=xzf^fl{r), dz=yzfz^,(r)^ 
d i= [(x 2 - y 2 )/2]/ 3 , 2 (r), and d 5 = [(3z 2 - r 2 ) /2V3]/ 3 , 2 (r) (see eqs. (125) - (129)), where 
d^ (C=l - 5) is a real function. We then have (d^\L\d^)—0 and (d^'\L\d^)= — (d^ \L\d{') 
for Ct^C; because L is purely imaginary and a Hermitian operator. Next, we con- 
sider ^j=J2c c ^d^, where ^ is chosen to be a real function. We have (^j\L\ i Sj) = 
Ec'Ec^M^C'I^K) = E C '>C C ? C ? (K'l-^K) + (d c \L\d c ))=0. Furthermore, when frf, 
we obtain 0M-^)=E C ' E^f cf (d c \L\d ( ) = £ c c ( p ^(d ( \L\d c )) =-^,\L\^ f ). 



54/55 



J. Phys. Soc. Jpn. Full Paper 

43) Regarding of eq. (79), we obtain C^^—C^^ using (\I>j/ |L M |\I>j) = — (tyj \L ^ ji) . A ' 2 In a similar 
way, 4°i= ^ is obtained using (<f 3 ,\L^ 3 ) = -{^ r ) . 42 From C'^=jC^ - £ { °l, we 

44) Using OMI^.) = -(^1^1^-0 (see ref. 42) and (^|^± |^) = (^|F^|^) of eq. 
(82), we obtain A±^= ^(l/ft) (*,|^|*i)/(^i - - -^ 2 (l/fi) 
E j( ^i)(*i|^l*j) - £y) = -A*^. In addition, we have (A±^)* = 
^ 2 (l/fi) E^d^iK^^I^-) (* i |(i A .)t|* 1 )/( J E 1 - = nX^l/h) E^d^iK"^)!^) 
(^•1^1*1)/^! - Ej) =-A^ = A% v , c 

45) J. Kondo: Prog. Thcor. Phys. 24 (1960) 161. 

46) See chap. 2 in ref. 31. 

47) In the case of a=c p = 0.15, the total sum of existing probabilities of d m and in i[>j of 
Eq. (122) is given by C 2 (£ 5 m=1{7tj) |c d | 2 + ^Li |c p | 2 )=0.14. Here, c p =0.15, which leads to 
C 2 En=i |c p | 2 =0.058, is roughly set by drawing upon the Fc 4p partial density of states (DOS) 
in the energy region of Fe 3d band of Fe4N. As to the DOS of Fe4N, see A. Sakuma: J. Phys. Soc. 
Jpn. 60 (1991) 2007. As to a schematic illustration of the partial DOS's of Fe4N, see S. Kokado, 
N. Fujima, K. Harigaya, H. Shimizu, and A. Sakuma: Phys. Rev. B 73 (2006) 172410. 

48) For example, see J. A. Stroscio and R. J. Celotta: Science 306 (2004) 242. In this literature, the 
typical vibration frequencies are considered to range from 0.01 THz to 10 THz. 

49) On the basis of the coefficient in eq. (7) in ref. 28, c q has been written as eq. (139). 

50) M. N. Leuenbergcr and D. Loss: Phys. Rev. B 61 (2000) 1286. 

51) A. Abragam and B. Blcaney: Electron Paramagnetic Resonance of Transition Ions (Clarendon 
Press, Oxford, 1970) Sect. 10.4. 

52) C. P. Slichter: Principles of Magnetic Resonance (Springer, New York, 1990) 3rd ed., Chap. 5. 

53) C. Kittcl: Introduction to Solid State Physics (John Wiley & Sons, New York, 1986) 6th ed., p. 
372. 

54) For example, see chap. 3 in ref. 53. 

55) For example, see M. T. Dove: Structure and Dynamics: An Atomic View of Materials (Oxford 
University Press, Oxford, 2003) Sects. 9.3 and N. 

56) For example, see J. M. Ziman: Electrons and Phonons (Oxford University Press, Oxford, 1960) 
Chap. I. 



55/55 



